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FLOOD DISTRIBUTION PROBLEMS BELOW OLD RIVER 


By Frederic M. Chatryl 


SYNOPSIS 


suring the safe passage of flood flows on the Mississippi River at and 
Old River requires that a controlled system of natural streams and ar- 
 floodways be operated to distribute peak flows in accordance with the 
1 capacities of the various system segments. Factors influencing the op- 
m include stages and discharges, sedimentation problems, economic and 
ercial considerations, and adaptation to constantly changing hydraulic 
‘ions. 


INTRODUCTION 


the limits of the watershed tributary to the Mississippi River are drawn 
ordinary map, the figure inscribed can be visualized as a large, upright 
, with its mouth stretching across much of the northern United States, 
3 spout in Louisiana. In this funnel the runoffis collected from a drain- 
isin 1,246,000 sq miles in area, in which are contained all or parts of 31 
states of the United States, and 2 Canadian provinces. Through the fun- 
spout pour 468,000,000 acre-ft of water in an average year, at rates 
ig from 100,000 cfs during low water season, to nearly 2,000,000 cfs in a 
1 high water, and to 3,000,000 cfs for the design flood. 

en the Flood Control Act of 1928 projected the federal government into 
31d of flood control in a preemptive manner, the mechanism which the 
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end of our funnel represents was relatively simple. Flow entered the 
via the Mississippi and Red Rivers, which streams are connected togeth 
a short length of channel called Old River, and was then conveyed to the 
of Mexico by way of the Mississippi River, and the Atchafalaya River. E: 
for levees, neither outlet was subject to any control, and flood manage 
was a matter of constructing and maintaining the levee system to suffi 
height and cross-section to insure its integrity. However, this was not e 
done, and levee overtoppings and failures were experienced in nearly ¢ 
significant high water. Before embarking on a description of the control 
tem that grew out of the Flood Control Act of 1928, some remarks on the 
eral topography and hydrography at and below the latitude of Old River a 
order. Present details of the area and of the flood control plan are shoy 
Fig. 1. Fig. 2 shows the area prior to start of construction of the flood 
trol plan. 


NATURAL SEGMENTS OF THE FLOOD CARRYING SYSTEM 
AT AND BELOW OLD RIVER 


The largest segment of the natural complex at and below Old River i: 
Mississippi. The Mississippi River adopted its present course below Old F 
approximately 1200 A.D., after abandoning the Teche and Lafourche dist 
taries. The present distance from Old River to the Gulf, via Southwest E 
is 321 miles. Below the gap at Old River, the Mississippi is fenced off t 
west by a continuous levee extending almost to the river’s mouth. On the 
bank, the river is confined by high bluffs as far south as Baton Rouge, ar 
a continuous levee thereafter. The carrying capacity of the leveed cha 
varies throughout its length. Just above Old River there is capacity 
2,700,000 cfs. This figure decreases to 2,100,000 cfs below Old River, tc 
proximately 1,550,000 cfs below Morganza Floodway, and finally to 1,25€ 
cfs below Bonnet Carre Spillway. The Mississippi River disposes of sli 
less than 70% of the total annual flow arriving at the latitude of Old River. 
though this percentage has been decreasing since 1855, exhaustive stu 
made in connection with the Old River control project do not reveal any 
nificant corresponding decrease in carrying capacity. As a matter of fact 
flow capacity of the River below Old River, as demonstrated by the st 
discharge relationship at that point, has been subject to no progressive ch 
at all, but has rather been subject to both increases and decreases in fl 
pacity. In 1950, for example, the stage-discharge relationship was as 
able as it has ever been throughout the period for which nip was 28 
available. , 

The second segment, with regard to size, is the Atchafalaya River. 
stream probably was spawned by excess water from the Mississippi 
spilling over the ridges of that stream’s abandoned Lafourche and 
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combined with overflow of cleared lands above the raft led to its removal 
dertaken in 1839. The removal was completed just prior to the War bet 
the States. 

Enlargement of the leveed portions of the Atchafalaya River during 
after the removal of the raft has been dramatically swift. At the time that 
moval of the raft was undertaken the normal discharge capacity at Simmes| 
La. (approximately 5 mi below the Atchafalaya head) was only 30,000 cfs 
1858, this figure had increased to 130,000 cfs. At the turn of the century 
figure had grown to 260,000 cfs and in 1950, the year of the last signifi 
flood at the latitude of Old River, it was 470,000 cfs. Even though there 
been no significant flood in the last 10 yr, the indicated normal capacity is 
approximately 500,000 cfs. It was this record of growth, in addition t 
fact that the Atchafalaya River gradient is double that of the Mississippi, v 
threw into sharp focus the imminence of what would have been an act of pi 
of monumental proportions, and led to the authorization of the $70,006 
project called Old River Control to prevent the abandonment by the Missis: 
River of its present channel below Old River in favor of the steeper rout 
the Atchafalaya River. 

Unfortunately for the American taxpayer in general, and flood control 
gineers in particular, the unleveed central reach of the Atchafalaya R 
which consists of a network of braided streams, many of which have a 
ments contrary to the general direction of the main flow, exhibited no ten 
cy toward self-enlargement. Instead, this portion of the river has con 
ously deteriorated through extensive siltation. The result of this siltatior 
been doubly detrimental to flow, not only reducing the available cross-sec 
but converting vast open water areas of low roughness to vegetated landa 
of high roughness. 

The Atchafalaya River terminates ina deep, well-defined outlet via Ber 
Bay and the Lower Atchafalaya River. Because of the small quantity o 
carried in this lower segment as a result of the sedimentation in the ce 
reach of the basin, this outlet is maintained in an efficient condition by th 
gress and egress of the tides. Eventually, however, the Atchafalaya wil 
gin to build a delta through its outlet section. In addition to this natural o1 
the Atchafalaya was provided an artificial outlet through Wax Lake in 194 

The smallest component of the overall system is Old River itself. 
stream, some 7 mi long, interconnects the Mississippi, the Red, anc 
Atchafalaya Rivers. Old River actually flows in one limb of an abanc 
Mississippi River bend. It is believed that sometime around 1500 A.D 
Mississippi River, in enlarging a westwardly migrating bend around wha 
later to become Turnbull Island, caved into the Red River, which was 
tributary to the Mississippi through what is now the Bayou Lettsworth <¢ 
nel. Subsequently, the Mississippi River continued its westerly meande 
with the result that Red River was forced to flow into the upper side o 
bend, At the same time, the junction probably resulted in diversion of 
away from Bayou Lettsworth and into a distributary near the present s: 
the city of Simmesport. This distributary is the present day Atcha 
River. 

After these conditions were established, the climactic event in the ia 
ment of Old River occurred. In 1831, Captain Henry Shreve, a pioneer i 
development of waterborne commerce in Louisiana, constructed the r 
which bears his name. Shreve’s cutoff severed the migrating bend previ 
referred to and formed Turnbull Island. The upper portion of the bend | . 


FLOOD DISTRIBUTION h) 


mown as Upper Old River, whereas the lower portion was called Lower 
ver. Shreve’s cutoff also halted the growth of the Atchafalaya River 
Jpper Old River no longer exists, except as a series of swales. Though 
elopment of Old River lagged that of the Atchafalaya for many years, it 
developing at a rate comparable to that of the Atchafalaya River. It is 
le that Old River is still inhibiting the growth of the Atchafalaya, inas- 
as it must continue to enlarge to meet the demands for water imposed 
Atchafalaya River at all stages. 

final segment of the natural system is the vast reservoir known as the 
ver Backwater Area. This is the area subject to inundation by Mississ- 
ver backwater through the gap in the levee system at Old River. Its 
embrace a total of 2,700 sq mi, and, at the elevation of the project de- 
90d, it is capable of storing 20,000,000 acre-ft of floodwaters. Because 
ubject to frequent overflow, the entire area is sparsely®8developed and 
inhabited. There is, however, a significant aggregate of improved lands 
1, some of which are now partially protected by levees, which makes the 
factor to be considered in distributing major floods. 


OVERALL SYSTEM FOR CONTROL OF FLOODS 
AT AND BELOW OLD RIVER 


ddition to the natural segments, the overall project includes several 
al segments that will be used in combination with the natural segments 
"ibute flood flows properly throughout the system. 

Atchafalaya Basin Floodway was formed by constructing “guide” levees 
d west of the Atchafalaya River, and generally parallel to the alignment 
stream. On the east side, the levee is carried upstream to approxi- 
the latitude of Melville, La., and then brought northeastward to join 
2 Mississippi River levee. On the west side, the levee extends above 
itude of Old River and ties to high ground in the vicinity of Mansura, 
ie floodway thus created is intended to provide one-half of the total ca- 
required to dispose of the project design flood. During such a flood the 
tte of inflow to the floodway would rise to as much as 1,500,000 cfs. In 
© introduce this total flow, three intakes will be employed. The first 
2 is the Atchafalaya River, which is not subject to control. The other 
2 artificial overbank floodways that flank the Atchafalaya River to the 
d west, both of which are subject to some measure of control. 
floodway to the east of the Atchafalaya River, known as Morganza 
ay, receives its flow directly from the Mississippi River through a 
ontrol structure. The structure is of reinforced concrete and consists _ 
sir 3,906 ft long equipped with vertical lift gates. The weir, with its 
t elevation 37.5 ft mean sea level is subdivided horizontally into 125 
ich of which provides a clear opening approximately 28 ft wide by 25 ft 
[The Morganza Floodway extends as far southward as the latitude of 
Jprings, La., where it will deliver its outflows to the floodplain of the 
Atchafalaya Basin. It averages approximately 5 mi in width. 

intake for Morganza Floodway is located on the Mississippi River 
7 miles below Old River, in the Raccourci-Old River area. Raccourci- 
rer is an ox-bow lake, formed by the Raccourci cutoff which occurred 
. The lake is expected to materially assist ingress of floodwaters into 
dway, and is further expected to act as a natural silt basin after the 
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The Morganza Floodway is designated to divert approximately 600,00 
from the Mississippi River to the Atchafalaya Basin. The lift gates pri 
complete flexibility in operation, in that flow through the floodway c:z 
started, stopped, or adjusted as required. The Morganza Floodway was | 
pleted in 1956, and has not been used since completion. An overall view c 
control structure is shown in Fig. 3. Details of the weir and gates can be 
in Fig. 4. 

On the western side of the Atchafalaya River lies the West Atchafa 
Floodway, larger in area than the Morganza Floodway and potentially n 
more controversial in operation. The West Atchafalaya Floodway differs { 
the Morganza Floodway in that control is not achieved by the positive me 
nism of gates, but by a levee of substandard grade called a fuseplug levee 
theory, the grade of such a levee is set to permit overtopping when the f 
in progress re&ches the point where operation of the floodway is require 
obtain the desired flow distribution. After overtopping, natural erosion o! 
levee, augmented as required by artificial degrading will be relied upon tc 
sure the required diversion. 

The West Atchafalaya Floodway extends southward to the latitude of K 
Springs where, like the Morganza Floodway, it will discharge its floodwa 
to the central Atchafalaya Basin Floodway. It averages approximately 8 n 
width. Unlike Morganza Floodway, which is practically uninhabited, the \ 
Atchafalaya Floodway supports a population of approximately 3,000 and « 
tains numerous improvements and much livestock. 

The West Atchafalaya Floodway has not been used, and will only come 
operation for a flood which approaches the project design value in magnit 
Fig. 5 shows the fuseplug levee at the head of the floodway. 

The final component of the Atchafalaya Basin Floodway is the Wax I 
Outlet into the Gulf of Mexico. This auxiliary outlet was constructed in | 
by dredging a channel from Six Mile Lake, through the Bayou Teche Ridge 
Atchafalaya Bay. The auxiliary outlet accomodates approximately 20% of 
total Atchafalaya Basin flow. An aerial view of the outlet is shown in Fig 

The third artificial floodway in the overall plan is Bonnet Carre Spill 
whichis located on the east bank of the Mississippi River 33 miles above | 
Orleans. It was constructed to divert excess Mississippi floodwaters 
Lake Pontchartrain. It is the oldest of the artificial floodways, having | 
completed in 1936. The control structure, which is located at the river en 
the Floodway, consists of a reinforced concrete weir, 7,698 ft in length, 
divided into 350 bays, 176 of which have a crest elevation of 15.8 ft msl, 
174 of which have a crest elevation of 17.8 ft msl. Control is by wooden t 
bers or “needles” which measure 8 in by 12 in. The Spillway was designe 
discharge a peak flow of 250,000 cfs at a stage of 20 ft at New Orleans. It 
been operated in 1937, 1945, and 1950. In 1945, the design discharge 
equaled or exceeded on 15 days. 

Significant problems have been encountered in the operation of Boi 
Carre Spillway, in that sedimentation in both the forebay area (the area 
tween the river bank and the control weir) and the Spillway proper are of 
ating to reduce its capacity in a progressive manner. Although the flow cz 
city can be maintained indefinitely by removing the silt deposits, such m: 
tenance is very expensive. View of the Spillway in operation are show1 
Figs. 7 and 8. 4 

Although it is not yet a functioning segment of the overall plan, the 
River Control project is sufficiently close to completion that consideratiol 
it should properly be included in any description of the overall system. 
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FIG. 5.—WEST ATCHAFALAYA FLOODWAY. SIMMESPORT-HAMBURG 
LEVEE AT HEAD OF THE FLOODWAY. 


FIG, 6.—WAX LAKE OUTLET. VIEWED IN THE DIRECTION OF FLOW. 


FLOOD DISTRIBUTION 2 


7.—BONNET CARRE SPILLWAY. VIEW DURING 1950 OPERATION, MISSISS- 
IPPI RIVER IN FOREGROUND. LAKE PONTCHARTRAIN IN BACKGROUND. 


: 8.—BONNET CARRE SPILLWAY. CLOSE-UP OF CONTROL STRUCTURE 
1950 OPERATION, ONE BAY PARTIALLY OPEN. 
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This project is one of the largest attempts at flood control engineering 
dertaken in the United States. Its basic purpose is to prevent the adoptio 
the Mississippi River of the shorter route to the sea offered bythe Atchafz 
River. The project includes a closure dike to permanently sever the na 
connection between the Mississippi and Atchafalaya Rivers via Old Riv 
gated control structure with inlet and outlet channels to provide for neces 
flow diversions from the Mississippi to the Atchafalaya under all condit 
and a navigation lock to preserve the navigational opportunities now affo 
by Old River. 

The control structure consists of two independent units in order tha 
necessary flexibility in operation may be assured. The so-called low 
portion consists of eleven openings, 44 ft wide. Of these eleven open 
eight are 55 ft deep and three are 70 ft deep. All openings are controlle 
multiple-leaf vertical lift gates, which will be manipulated to control lowy 
and medium stage flow. The overbank or high sill structure consisting | 
bays, 44 ft wide and approximately 12 ft deep, equipped with needle type 
sures, will control diversion of high water and major flood flows. Ths 
structures are capable of passing a flow of approximately 700,000 cfs u 
existing channel conditions in the Mississippi and Atchafalaya Rivers. | 
tinued development of the Atchafalaya River will increase this maximum 
pacity. 

At the present time, the control structure, the outflow channel, and mo 
the levees are complete. The navigation lock is under construction, an 
entire project is scheduled for completion in 1964. However, should an er 
gency develop in which the need for closing Old River became immediate 
control structure could be operated. 

In general, Old River Control will be operated in such a fashion as to 
tribute the flows artificially in almost the same proportions in which 
would have been distributed under the natural conditions that existed in | 
The details of operation are not appropriate for presentation herein; howé 
the entire operation will be oriented around the proposition of preservin 
existing stability of the Mississippi River channel below Old River. An c 
all aerial view of the salient project features is shown in Fig. 9. 


OPERATION OF THE VARIOUS FLOODWAYS FOR DISTRIBUTING? 
FLOOD FLOWS AT AND BELOW OLD RIVER 


The natural elements of the flood control plan below Old River are | 
operating. The artificial floodways and control structures are operate 
sucha manner as to insure that the safe capacities of the various segmen 
the system are not exceeded. In the operation, guide lines are provide 
stages and discharges. There is only one specific operational criterion 
the sanction of law, and it is that the Bonnet Carre Spillway shall be ope! 
to prevent stages at New Orleans from exceeding 20 ft. 

The natural leveed channel of the Mississippi River at and below | 
Orleans can accommodate all floods in which the peak flow in that reach c 
river does not exceed 1,250,000 cfs. A flood producing this peak flow c: 
expected to occur on an average of once every 5 to 7 yr, and would produ 
maximum stage at New Orleans of 20 ft. Such a flood represents the a 
flood insofar as floodway operation is concerned, 


FLOOD DISTRIBUTION 11 


r floods exceeding the magnitude of that previously described, operation 
» or more of the artificial floodways is required. In a typical flood the 
ice of operations would generally be as follows. 


Operation of either Bonnet Carre Spillway or Morganza Floodway when 
in excess of 1,250,000 cfs is in sight for the Mississippi River below 
iver. 

Operation of the remaining Mississippi River floodway when a flow in 
3 of approximately 1,500,000 cfs is in sight for the Mississippi River 
Old River. 

Operation of the West Atchafalaya Floodway when the total flow at the 
le of Old River exceeds approximately 2,800,000 cfs. 


OPERATIONAL PROBLEMS 


he relationships between stage and discharge throughout the various seg- 
of the flood control system below Old River were constant, the opera- 
procedure would be straightforward and invariable. Whereas flow in 


).—OLD RIVER CONTROL. LOW-SILL CONTROL STRUCTURE AND OUTLET * 
_ CHANNEL IN CENTER OF VIEW. HIGH-SILL CONTROL STRUCTURE TO 
_ THE RIGHT. | 


l rivers follows the same immutable laws governing all fluid flow, 
al changes in the hydrography of these streams produce continuing 
s in their hydraulic efficiency. In addition to well-defined changes of a 
ssive nature, such as those taking place on the leveed portion of the 
Jaya River, other changes ‘for which cause-and-effect relationships 

been established combine with the aforementioned, and operate to 
n inflexible operational procedure impracticable. Notable among these 

iation in stage-discharge relationships for the Mississippi River it- 
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The intake capacities of all the artificial floodways are affected by 
changes previously mentioned, The Morganza Floodway intake capaci 
indicated to be in the neighborhood of 50,000 cfs per ft of stage for pr 
design conditions. Inasmuch as the stage-discharge relationship for the 
sissippi River below Old River has varied by as much as 2 ft for large f! 
the extent to which these variations influence diversions through Morg 
Floodway can be readily appreciated. The significance of these variatio 
enhanced when it is considered that they tend to vary from one section 0 
river to another. For example, the stage discharge relationship in the 
sissippi River between Old River and Baton Rouge was, during the 1950 f 
as low as it had been for the entire period of record, whereas such wa: 
the case below Baton Rouge. 

The performance of the West Atchafalaya Floodway is influenced by ) 
ations in stage-discharge relationships to an even greater extent. For e: 
ple, as recently as 1941, the planned distribution of the project design 
anticipated simultaneous flows of 650,000 cfs and 250,000 cfs, inthe Atchaf 
River and the West Atchafalaya Floodway, respectively. In 1950 the : 
along the fuseplug levee at the head of the floodway averaged 5.5 ft belo 
fuseplug levee grade eventhough the peak flow inthe Atchafalaya River rez 
650,000 cfs. The intake difficulty imposed by the steady lowering of the st 
discharge relationship for the Atchafalaya River at Simmesport, La. w: 
worsened as time goes by. 

The general problem of properly distributing floods below Old Rive: 
volves providing answers to the following questions. 


1, At what time during the flood should the first floodway be oper ated 

2, Which floodway should be operated first? 

3. At what time during the flood should the second floodway be oper 

4. If the flood is indicated to be of extreme magnitude, what steps s 
be taken to insure that the West Atchafalaya Floodway operates in a ti 
and effective manner? | 


Time of Initial Operation.—This problem involves careful considerati 
antecedent stages and discharges on the Mississippi River. Because pr 
stages will be approached along the entire Mississippi below Bonnet ( 
Spillway prior to the operation of the first floodway, whereas Atchaf 
Basin stages will be relatively moderate at that time, control and limit 
of stages in that reach or river are the primary concern. Accordingh 
operation must be closely coordinated with the stage at New Orleans. | 
narily, a fairly definite decision can be reached as to the proper timing ¢ 
initial operation. Although unusual, cases will arise in which the proper 
of initial operation will be somewhat obscure, This could happen in a 
cresting at a discharge slightly larger than 1,250,000 cfs in the Missis 
River, particularly if the flood hydrograph rose rapidly. Such a flood 
probably produce unusually low stage-discharge relationships on the Mi 
sippi River in the early phases of the flood, in which case there woul 
tendency to delay the operation in the hope of avoiding it altogether. This 
cedure could result in waiting too long, and thereby permitting the New Oi 
stage to go above 20 ft. On the other hand, early operation based on 
charges alone, might result in an operation that could otherwise be av 

No all-inclusive specific rule canbe laid down for time of initial a 
but, in general, it would appear proper to base a Morganza first operati 
Mississippi River discharge, and a Bonne Carre first operationon New Or 
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_ Thus, if Morganza Floodway were to be operated first, this operation 
come sufficiently in advance of the Mississippi River discharge reaching 
000 cfs to insure that not more than 1,200,000 cfs would flow past New 
ns. If Bonnet Carre Spillway were to be operated first, the operation 
come at a stage of approximately 19+ ft at New Orelans, with more than 
n prospect. 

ich Floodway to Operate First.—This question is agreat deal more com- 
han that of when to operate. There are advantages and disadvantages to 
perational sequences, The operation of Bonnet Carre Spillway is physi- 
quite simple. The land on which the spillway is located is owned in fee 
= by the federal government, and the problem of evacuation is nil. The 
t Carre operation involves little or no economic or commercial dis- 
n. Although unanimity of opinion does not exist with respect to the ef- 
fa spillway operation on fish and wildlife opportunities, most informed 
‘S agree that operation of the spillway is helpful in this regard. Based 
erienced flows in Bonnet Carre Spillway and theoretical computations 
iodel studies for Morganza Floodway, it is indicated that the magnitude 
flood that the overall system can safely accommodate with only Bonnet 
Spillway in operation is greater thanthat which would be accommodated 
mly Morganza Floodway. In practical terms, this means that floods 
ngnear 1,500,000 cfs inthe Mississippi below Old River could be handled 
sonnet Carre alone without exceeding 20 ft at New Orleans, but not with 
nza alone. 

2 most serious disadvantage of a Bonnet Carre operation lies in the sed- 
ition concommitant to operation. In addition, the stage lowerings real- 
‘om a Bonnet Carre operation are limited essentially to the Mississippi 
below the Spillway itself. 

» Morganza Floodway, on the other hand, provides stage lowerings and 
‘S in the Red River backwater area as well. Based on present best esti- 
_the sedimentation problem in Morganza Floodway will be much less 
s than that in Bonnet Carre Spillway. 

the unfavorable side, operation of the Morganza Floodway is much more 
2x than that of Bonnet Carre Spillway, and its effects will be felt over a 
r area, There is a substantial evacuation problem involved. Backwater 
g within the West Atchafalaya Floodway Area will be increased. Drain- 
the Pointe Coupee Protected Area, and in the towns of Melville and 
Springs, will be interrupted or impaired, and although dedicated areas 
ding runoff have been provided, there is always the outside chance that 
‘eme rainfall will overtax the storage provided. 

an overall basis, it would appear that the advantage lies with a Bonnet 
first operation. However, it must be borne in mind that sedimentation 
211 be a controlling element. In any event, it is planned to operate the 
za Floodway first in the next flood requiring operation of a floodway, 
r that the definite information needed tocalibrate the control structure, 
ine overall hydraulic effects, and evaluate the sedimentation problem 
antitative manner may be collected. 

e of Operation of Second Floodway.—The factors governing this prob- 
e similar to those relating to time of operation of the first, except that 
eshold discharge will be approximately 1,500,000 cfs. As before, if 
za is to be operated, the operation would be keyed generally to the 
‘ge, whereas a Bonnet Carre operation would follow more closely the 
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Operation of West Atchafalaya Floodway.—Although the West Atchaf 
Floodway is not subject to control in the same sense that Morganza Floc 
and Bonnet Carre Spillway are, it nonetheless presents a most comple 
erational problem. The central element in this problem is the leveed ch 
of the Atchafalaya River. As mentioned previously, the operation of the 
Floodway is automatic and will result from the stage along its fuseplug | 
exceeding the levee grade. The natural enlargement of the Atchafalaya 1 
has progressively lowered the stage along the fuseplug for a given total 
at the latitude of Old River. Although such a development may be view 
generally favorable, in that it is producing a continuous increase in the « 
all intake capacity at the head of the Atchafalaya Basin, it does introdu 
impediment to timely operation of the West Atchafalaya Floodway, be: 
such operation must await the natural overtopping of the fuseplug. Failv 
the West Atchafalaya Floodway to operate on time would probably res 
overtaxing the safe leveed capacity of the Mississippi River, particular], 
low New Orleans. 

At the present time, it appears that operation of the West Atchaf 
Floodway will commence early enough in a great flood to limit Missis 
River flows to a safe value, provided that the initial overtopping of the 
plug is followed promptly by various measures designed to facilitate ent 
floodwaters into the Floodway. Such measures might include degrading | 
fuseplug levee at the head of the floodway or its destruction by explos 
Similar action may be taken along the West Atchafalaya River levee, i 
quired. However breaching the West Atchafalaya River Levee shou 
avoided, if possible, because it would produce excessive velocities i 
Atchafalaya River above the breach, and could result in uncontrolled 
caving and loss of the bridge across the Atchafalaya River at Simme: 

In view of the fact that continued improvement of the Atchafalaya Ri 
. expected, the West Atchafalaya Floodway problem will tend to grow pro 
sively worse. The long-term solution lies in revision of the grades of 1 
protecting the floodway. 


CONCLUSIONS 


Any plan for the control of floods below Old River must possess the fe 
ing attributes: 


1. It must be capable of safely distributing the maximum peak flow: 
may reasonably be expected to occur. 

2, It must accomplish this distribution with a minimum of economi 
ruption, 

3. It must be adaptable to changing conditions, so that its effectiv 
may remain unimpaired in the future. 


The system described herein is shown to satisfy these requiremen 
can be operated in a planned and predictable manner for peak total fl 
the latitude of Old River up to and including 3,030,000 cfs, a value 
been established by extensive meteorological iid hydrological audi 
largest likely to occur. With the exception of the West Atchafalaya Floc 
which will operate less often than once every hundred years, the economi 
ruption deriving from system operation is not great. It varies from neg] 
in the case of an operation of Bonnet eine Spillway to minor for a Mor 
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y Operation. And the overall system’s flexibility and adaptability are 
t the end of its useful life cannot now be foreseen. With the system 
the threat of damaging overflow below Old River has been virtually 
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DICTING STORM RUNOFF ON SMALL EXPERIMENTAL WATERSHEDS 


By Neal E. Minshall,! M. ASCE 


SYNOPSIS 


reful estimates of rates and amounts of runoff are required for design 
od retarding and channel stabilizing structures, and bridges and culverts 
all upstream watersheds. The period of record of many of the small ex- 
1ental watersheds in the United States is too short topermit direct analy- 
r reliable estimates of the magnitude and frequency of storm runoff. A 
dis presented for extending the period of runoff records based on analy- 
existing short term records of rainfall and runoff for the watershed, and 
ser record of the rainfall alone. The method involves: (1) estimating 
| runoff volumes from the rainfall pattern and antecedent rainfall, and 
stributing this runoff through an adaptationof the unit hydrograph princi- 
\ method is also presented for developing synthetic unit hydrographs for 
ed areas. 


INTRODUCTION 


yery limited amount of runoff data is available on small watersheds, and 
records are not generally of sufficient length to provide reliable esti- 
for design floods. Since records of precipitation are generally avail- 
ar much longer periods than runoff records, it would be desirable to es- 
n precipitation-runoff relations and use such relations to develop flood 


\ 


2,—Discussion open until January 1, 1961. To extend the closing date one month, 
en request must be filed with the Executive Secretary, ASCE, This paper is part 
opyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
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hydrographs from the longer precipitation record. The method presented 
is summarized as follows: 


(1) Develop time-retention curves for each soil cover combination, bh 
on an antecedent precipitation index and season. 

(2) Estimate the storm runoff volume from the precipitation pattern an 
time-retention curve. 

(3) Derive unit hydrographs for various rainfall intensities. 

(4) Apply these relationships to a number of outstanding rainfalls and « 
pare the computed hydrographs with the observed runoff. 


The relationships developed can be applied to a longer precipitation re 
from a nearby recording rain-gage. A relationship has also been develop: 
permit adaptation of these principles to ungaged areas. 

The analysis for developing the procedures was made withrainfall and 
off data from three small experimental watersheds located near Edwardsy\ 
Ill. (Fig. 1). 


DESCRIPTION OF EXPERIMENTAL WATERSHEDS 


The three Edwardsville watersheds, with a large number of runoff pei 
each year, offer an excellent opportunity for studies of rainfall-runoff 1 
tions on small watersheds. These watersheds (details shown in Fig. 2) 
designated as: W-1, a 27.2-acre cultivated area; W-2, a 50-acre past 
area; and W-4, a 290-acre mixed cover watershed. The two small waters 
are part of the 290-acre area. These watersheds are, in general, quite 
in the upper portions but drop off rather abruptly near the waterways. 

Watershed W-1 is a fan-shaped area with a range in elevation of 20 ft, 
more than two-thirds of the area has slopes less than 1%. This area was 
cultivated during the years of 1938 and 1948-50, 80% cultivated in 1951 
1955, 100% alfalfa in 1940 to 1943, and 44% cultivated during the rest o: 
period of record. Areas not cultivated were in legume hay. 

Watershed W-2 is an oval-shaped area having a length about one and 
half times its width. The difference in elevation between the divide ant 
runoff station is 39 ft. About one-third of the area has an average slo} 
1.5%, and another third has 12% slopes. During the period of record, cove 
this area has varied from 0 to 15% cultivated with the remainder in pastu 

Watershed W-4 is nearly circular and has a range inelevationof 63 ft. | 
third has anaverage slope of 1%, and another third has 12% slopes. The n 
mum slope near the waterway is 35%. There are four small ponds with a 
drainage area of approximately 30 acres or 10% of the watershed area. — 
ing the period of record, from 17 to 30% of this area has been cultivates 
the remainder mostly in pasture. 

The soils on these watersheds are of loessial origin over glacial till 
an impervious claypan layer at depths of 10 to 20 in. These soils, whel 
and protected with a good vegetative cover, will take in precipitation r: 
rapidly until the surface layer becomes saturated, after which additional 
cipitation willresult in a high percentage of runoff. The data obtainedon { 
areas are applicable to the central claypan areas of southern Illinois, s 
Indiana, southwestern Ohio, Missouri, eastern Kansas, and northeastern ¢ 
homa. The mean annual precipitation over most of this area is about 4 
but the mean monthly rainfall varies widely from a fairly uniform distri 
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FIG. 1.—LOCATION OF WATERSHEDS 


20 August, 1960 ] 
in southern Illinois, to eastern Kansas, which has high spring and summe! 


low winter precipitation. 


BASIC DATA 


Continuous records of precipitation, runoff, temperature, humidity, and 
cover were collectedon the three Edwardsville watersheds from March, 1 
through December, 1955. 


—-—~- Waterway ~= 


— — — Watershed boundary nt / 
—500— Contour mae 7 


&, Concrete drop spillway F 
Cee Farm pond 0) 400 800 1 
@ §Gaging station . : 

e Rain gage location 1938-1955 Scale in feet 7 
Rain gage location 1938-1947 Contour interval 10 feet 
fe} Rain gage location 1947-1955 


: 
FIG, 2,—MAP OF EDWARDSVILLE, ILL., WATERSHEDS 
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recipitation was measured with seven recording raingages from 1938 to 
and five recording gages from 1948 through 1955. Additional recording 
$ were in operationon the two small areas during anintensive infiltration 
-from 1940 to 1943. The only gage on the 27-acre area, except for 1940- 
vas located near the runoff station. Because of the variation in rainfall 
ig thundershowers, this gage may not be representative of the actual aver- 
recipitation on the area. 

inoff from these areas was measured with concrete V-notch broad crest- 
irs equipped with water stage recorders. Runoff records from watershed 
are based on the effective area varying in size from 290 acres, when all 
5 are at spillway elevation prior to rainfall, to 260 acres when noneof the 
3 are discharging through the spillway. 

ofinitions.— The terms used herein are defined as follows: 


Retention - That part of the precipitation not appearing as runoff, in- 
S plant interception, depression storage, and infiltration. 

Antecedent Precipitation Index (API) - An empirical measure of the ef- 
Mf precipitation falling a given number of days prior to the period of ex- 


Runoff - As used here includes all flow that eventually passes over the 
including surface runoff, interflow, and base flow. 

Unit Hydrograph - A hydrograph representing 1 in: of runoff. 

Time to Peak - The time from beginning of excess rainfall to the peak 
9 runoff for storms of short duration. 

Period of Excess - Thetime during which precipitation occurred at rates 
sess of the infiltration capacity. 

Infiltration Capacity - The maximum rate at which a given soil, when in 
2 condition, can absorb falling rain.2 


ESTIMATING STORM RUNOFF VOLUMES 


alysis was first directed toward developing a relationship between ante- 
t precipitation (as a measure of the available soil moisture storage) and 
tion. 

mputing Antecedent Precipitation Index.—The length of time that must be 
Jered in computing any antecedent precipitation condition will vary with 
‘rmeability of the soils and subsoils, temperature, evaporation, transpi- 
, and vegetation. For very permeable soils such as sands or the better 
ams, a 5-day antecedent precipitation index may suffice, whereas, for a 
y permeable soil it may benecessary touse a periodof 30 days or longer. 

us methods and lengths of period were tested in computing the antecedent 
vitation index for Edwardsville areas, and the one that gave the best corre- 
with retention was a 30-day index plus 20% of the departure from normal 
30- 60 day precipitation. The final selection which gave nearly as good 
on with less computation was basedon a 30-day periodusing the formula: 


API = Po + Pyk + Pok2 + Pgk9 + - --- - Pyk”........ (1) 


ch API is the antecedent precipitation index, Py refers to precipitation 
24 hr prior to the excessive period which produced the runoff shown in 
SS a; eee ae 
Surface Runoff Phenomena, Part I—Analysis of the Hydrograph,” by Robert E, Hor- 
iblication 101, Horton Hydrological Lab., Voorheesville, N, Y. 
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column 6 at Table 1, Pj, Po, and P3 indicate precipitation 1, 2, 3 days, anc 
on, prior to the storm, n denotes the number of days used to establish the 
dex (in this case 30), and k is a constant depending primarily on the soil t 
(0.95 gave the best results for the claypan areas). 

Relation of Antecedent Precipitation to Retention.—Initial investigati 
showedthat a large volume of data was needed to develop relationships betw 
antecedent precipitation index andretention. Therefore, consideration was gi 
to combining the data from watersheds W-1, W-2, and W-4. Comparisol 
amounts of precipitation retained from the excessive portionof storms sho\ 
that retention on W-1 (27.2 acres - all alfalfa), during 1940-1943 was ak 
95% of that on W-4. When W-1 was 96% cultivated, the retention for individ 
storms averaged 18% less than W-4. Retention on W-2 (50 acres - of all p 
ture) averaged 5% more than on W-4. Comparison of annual and seasonal r 
offs on watersheds W-1 and W-2 against W-4 show that for the period of 19: 
1955, the average runoff on the two small areas was only 9% more than on 
4. With the apparent small variation between watersheds, it was felt that d 
from all three areas could safely be combined for the retention analysis. 

Table 1 givesthe data for watershed W-4 for allstorms during the 18 y1 
record that caused runoff in excess of 0.1 in. and where the precipitation 1 
such that all runoff could be charged to a definite period of rainfall exce 
Where it was possible tobreak a storm down into several distinct runoff pe 
ods, each of these periods was entered in the table. Some winter storms w 
eliminated because of frozen ground or melting snow. Similar informatior 
that presented in Table 1 was. compiled for watersheds W-1 and W-2. 

The data for all watersheds were divided into several groups according 
season and further subdivided by groupings basedon the API. In grouping th 
data according to season, one should recognize that variation in season fr 
one year to the next may be as much as two weeks. Because of this, it is 
entirely logical to arbitrarily select a given date as the limit of a particu 
group. The year was first divided into spring, summer, fall, and winter v 
three subdivisions for API in each, giving a total of 12 groups of data. 

The data from columns 4 and 6of Table 1 andsimilar information from ¥ 
and W-2, grouped as explained previously, were plotted on rectangular coor 
nate paper. The results indicated a rapidly changing infiltration capacity d 
ing the spring and fall seasons and fairly constant conditions during the su 
mer and the dormant season. They also indicated the advisability of slig] 
different grouping according to seasons. The final groupings, selected as ¢ 
ing good results with a minimum of groups, were: November 15 to April 
April 15 to May 15; May 15 to June 30; July through September; and October 
to November 15. é 

Fig. 3(a)shows these time retention curves for the period July through s 
tember at different antecedent precipitation indexes. Curves were fitted i 
the groups of points as plotted. Similar sets of curves were prepared for e 
of the periods listed. Fig. 3(b) shows the maximum retentions that result fr 
antecedent precipitation indexes less than 0.5. Fig. 3(c) shows retentions 
the various periods, termed average because they occur more frequently 
the others. To continue these curves to the origin of time and depth z 
serve no purpose. Those representing dry or summer conditions would br 
down very sharply, indicating the magnitude of plant interception, depress 
storage and initial infiltration rate. Minimum infiltration rates on these st 
are nearly constant throughout the year varying between .02 and .06 in. per 


Date 


9, 1939 


, 1940 


7, 1941 


1942 
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TABLE 1.—DATA FROM SELECTED STORMS CAUSING 
MORE THAN 0,1 IN, RUNOFF, 1938-55 


Precip. Excessive Runoff Retention 
before Precipitation from from 


excess Amount,/ Duration excess, excess, 
Inches Inches |Hrs. Min.] Inches Inches 


.04 85 5 30 19 .68 
0 .76 2 40 44 32 
0 .59 2 07 =o 31 
0 45 0 28 39 06 
0 .66 0 22 .20 46 
0 158 3 50 as 40 
10 68 1 45 152 16 
.04 1.30 0 32 63 67 
0 .71 0 08 .20 51 
0 68 0 55 16 52 
93 3.28 2 36 1,33 1.95 
53 36 0 13 at 09 
.03 93 5 30 ay 56 
15 61 2 17 18 43 
0 .40 1 50 127 13 
0 155 4 00 33 22 
0 43 5 00 13 30 
55 1.98 1 20 .58 1.40 
08 .84 0 12 22 62 
0 .58 0 40 19 39 
25 AT 0 45 16 31 
07 4,42 9 20 1.80 2.62 
05 1.14 4 07 135 79 
05 3.58 2 35 “73 2.85 
53 A9 0 13 .16 33 
0 1.25 0 55 (25 1.00 
0 1.43 1 00 ee 1,22 
0 .57 0 25 .24 133 
.07 64 0 17 .10 .54 
.04 49 0 20 12 37 
.59 54 0 40 .29 125 
0 1.95.98) 17 00 .78 1,17 
10 1.00 1 45 19 81 
.05 1.16 0 33 .40 .76 
0 1.39 0 25 .69 .70 
05 1.05 0 20 7 134 
05 1.70 4 15 .99 71 
£05 2.50 1 20 1.78 wr 
.05 3.85 7 25 2.37 1,48 
2.50 .75 0 30 152 .23 
0 122 1 50 .09 13 
0 er 1 30 .57 .24 
0 .56 2 25 12 44 
0 1.06 0 23 _.33 73 
0 ay 0 25 ail .46 
.07 .69 0 18 125 44 
.20 2.31 6 50 | 1.64 67 
0 .60 7 apenas Bl: .43 
0 48 3 10 og 125 
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Precipita- 
tion Index 
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TABLE 1 (CONT’D) 

(1) (2) (3) (4) | (5) (6) (7) 
Mar. 16, 48 134 1 20 124 10 1.54 
Mar. 19, 0 52 3 30 .36 16 1.85 
May 6, 95 1.00 2 10 .39 61 1.9% 
May 8, 0 £30 i 30 .12 .18 2.8( 
May 10, 0 41 1 35 .36 .05 3.61 
May 15, 35 135 0 12 15 .20 3.91 
May 16, 0 .29 0 07 16 .13 4.03 
May 17, 0 5.05 | 29 00 3.96 1,09 4,3( 
May 17, .50 4,55 | 24 00 3.55 1.00 4,7 
May 17, .50 1.85 6 40 1.44 41 4,7 
May 17, 2.65 1.00 1 45 .81 .19 6.5 
May 17, 3.70 1.65 7 40 1.29 36 7.51 
May 19, 0 .60 8 15 AT .13 8,1! 
June 5, 0 .87 2 45 ia .75 4A 
June 6, 0 1.43 3 30 .82 61 5.0. 
June 10, 0 std 0 50 .34 .o9 4.8 
Apr. 10, 1944 .06 .53 1 05 .08 45 11 
Apr. 22, 509 4.67 |14 00 3.15 1.51 2.1) 
Apr. 22, 0 1.01 4 00 .68 .33 6.9 
May 22, 0 1.06 2 10 19 .87 2.4 
May 23, 0 .59 1 10 .13 .46 2.3. 
May 23, 0 41 1 40 .13 .28 2.9 
Aug. 30, 0 1.70 2 00 “22 1.48 1.9 
Feb,21, 1945 Bie: 185 4 45 .18 .67 8 
Mar, 1, 0 1.49 |16 20 59 .90 1.8 
Mar, 6, 0 .64 4 00 .38 .26 3.9 
Mar, 25, 0 1.78 ob 18 30 .75 1.03 2.1 
Mar. 29, 0 1.32 1 55 61 71 2.9 
Mar, 30, 0 1.43 21 00 .96 47 4.3 
Apr. 1, 0 54 6 15 134 .20 5.5 
Apr. 12, .10 1.31 5 00 42 .89 34 
June 9, 0 52 0 13 Rit .25 3.1 
June 12, 0 Bz 1 30 .33 .54 3.2 
June 16, 0 2.20 1 50 1.33 .87 3.1 
June 16, 0 So 0 18 .18 14 5.4 
June 17, 0 124 2 00 .14 .10 5.6 
June 17, 0 134 4 40 .20 114 5.8 
Oct earl; 0 1,44 6 30 ae 1,22 1.0 
Nov, 12, .05 .64 0 40 .16 48 1.0 
Nov. 26, 0 .57 1 35 .12 45 ld 
Dec, 24, 0 oe gd ht eg 30 33 Bi 3 
Jan, 11, 1946 ) Ge 4 45 24 29 isa 
Feb, 12, 0 Tee 21 00 87 85 4 
Mar, 5, .10 1.40 5 00 74 .66 14 
Mar. 17, 0 aH 14 00 65 Ae ae 
May 3, 0 .64 0 30 “27 iBT 2.2 
May 17, 0 1.05 6 05 35 70 28 
Aug. 5, 0 1,12 0 38 17 .95 2.4 
Aug, 14, 15 1.49 3 40 132 ale hd 4,0 
Aug, 14, 20 3.09 2 20 2.31 .78 Bp 
Aug, 14, 20 5.00 6 05 3.41 1,59 2 
Aug, 14, .20 1.10 | 15 .67 43 5. 
Aug. 15, 0 1.45 0 50 1.10 .35 9, 
Aug. 15, .30 3.66 | 10 40 2.54 12 10. 
Oetaai, 5) 1.49 6 30 114 1,35 “1 
Oct. 31, 0 4.01 S15 45 | 1.75 2.26 2. 
Nov. 2, 0 91 4 20 .50 Al 6 
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1948 


1949 


1950 


1951 


25 
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TABLE 1 (CONT’D) 

(2) (3) (5) (6) (7) 
0 1.05 9 30 44 61 6.40 
0 1.45 8 00 .89 56 6.40 
.06 .50 5 00 15 35 5.59 
0 1.00 6 40 70 .30 1,31 
0 .93 3 25 AT 46 Ao 
0 48 4 30 25 23 py 
0 53 3 30 .39 14 1.57 
0 81 7 05 .33 48 1.70 
0 1.59 15 00 98 61 LRT 
0 1:37; co 43 95 42 3.16 
0 55 3 55 10 45 sis 
0 91 6 00 34 57 96 
0 1s ll 20 59 54 1.78 
0 57 1 35 30 Dr 2.39 
0 .70 6 50 .30 .40 2.68 
0 26 2 00 7 .09 3.20 
0 .76 1 20 18 58 2.04 
0 84 0 50 au, ie, 1.76 
0 59 0 30 E83 46 2.63 
0 54 0 10 10 44 3.07 
) 52 0 30 14 38 3.55 
0 85 1 10 35 50 4,23 
0 2.45 6 10 1.10 15) 4,59 
0 1.85 1 45 85 1.00 4,59 
0 .98 4 05 34 64 02 
0 1.81 12 00 1.01 .80 1.34 
0 81 0 15 13 68 .98 
0 1.10 0 23 58 52 1.78 
0 2.02 2 10 68 1,34 230 
15 1,35 1 05 1217 1.08 1,18 
04 1,00 0 30 Al 59 2.30 
0 1.86 1 10 49 1 S7 1.47 
45 1.67 4 00 22 1.45 1.83 

1.10 85 5 00 99 63 3.85 
0 42 0 20 .10 32 2.89 
0 59 4 15 25 34 2.70 
0 4,66 25 00 3.65 1.01 2.23 
0 59 5 30 46 13 5.19 
0 55 2 40 13 42 2.89 

25 1.07 5 10 .60 AT 1.30 
0 54 5 00 33 aor 2.33 
0 52 6 40 32 20 2A 
.10 74 2 50 19 55 1.54 * 
0 86 at 20 24 62 1.57 
0 1.50 6 30 62 88 2.35 

50 .98 3 00 51 47 2.85 
0 127 1 30 .33 94 2.49 
0 62 1 40 26 36 3.58 
0 52 3 00 .19 33 12 
0 1.65 15 30 70 95 Pe 
0 174 2 00 62 er 2.42 
0 .85 1 45 14 Brg 3,88 

.10 50 0 25 ait 39 4,55 
0 92 0 13 2 ht 4,42 
0 92 0 35 43 .49 4,46 
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TABLE 1 (CONT’D) 
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Mar, 3, 1953 


DMhrowhH bw 


Feb, 19, 1955 


0 
0 
0 
0 
0 
5 
0 
2 
0 
3 
" 
1 
2 
3 
0 
0 
1 
0 
0 
2 
7 
4 
5 
0 
i 
2 
0 
1 
0 
0 
5 
0 
2 
3 
2 
0 
0 


These studies show that differences in antecedent precipitation have ct 
siderable influence on retention during the summer but are of minor importa 
in the winter. The best correlation between antecedent precipitation index < 
retention was obtained for the period of October 15 to November 15, but t 
might change markedly if a larger number of points were available. The 
curves are based on the assumption that, for a particular season, the infilt 
tion capacity is a function of antecedent soil moisture and is independeni 
rainfall intensity. Because these are mass retention curves, the slope o 
tangent is a measure of the infiltration capacity rate. 

Because it is the amount of precipitation retained rather than the total p 
cipitation that adds to the soil moisture, a refinement was attempted usin 
antecedent retention index. The results showed no noticeable improveme 
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The investigations also indicate that data from a 187-acre watershed 
Hamilton, Ohio, and the 153-acre area on the McCredie, Mo., experiment 
tion agree quite favorably with the retention curves developed fromthe Edwa1 
ville data. The soils on these areas areof thesame general type as at Edwa 
ville, but depths to the impervious claypan may be somewhat different. 

Testing the Retention Curves.—Totest the accuracy of the retention cur 
the estimatedrunoff volumes for a number of storms were obtained by pla 
the time retention curve for the appropriate season and antecedent precij 
tion index under the mass rainfall diagram, with the zero time and deptha 
beginning of rainfall excess. In the case of watershed W-4, the mass rai 
diagram was the average of all gages. The vertical distance between the 
retentionand mass rainfall curves represents the estimated amount of rai 
excess, up to any time, which will eventually pass over the weir. Propel 
lowance must be made for the decrease in depression storage and plant ir 
ception, when, during any periods within the storm, the rainfall ceases or 
below infiltration capacity. 

A large number of selected storms that gave an estimated runoff of 1/ 
or more were tabulated and compared with the observed runoff. These ¢ 
parisons for watersheds W-1 and W-4 are shown in Figs. 4(a) and 4(b), 
spectively. In general, the best agreement was obtained for storms occur 
during the dormant season and at other timesof high antecedent moisture 
ditions, such as July 9, 1942, May 17, 1943, and August 15, 1946. Alth 
individual values show considerable variation between the estimated and a 
runoff, the totals for all stormson W-1 agreevery well. The estimated amc 
on W-4 appear to be generally too high, indicating the possible need for ad: 
ment of the time retention curves on this area. 


DEVELOPING HYDROGRAPHS 


After establishing a method for estimating storm runoff volume, the 
step in the analysis was todevelop a method of estimating the shape of the 
drograph. 

Rainfall Intensity and the Unit Hydrograph.— Attempts to establish an a 
age unit hydrograph for each of the Edwardsville watersheds showed a 
variation between storms having different rainfall intensities. It was desir. 
therefore, to establish relationships between rainfall intensity and shape o 
unit hydrograph. = 

For the 27-acre area, a seriesof storms was selected having isolated : 
ods of intense rainfall of 10 to 15 min duration, of reasonably uniform inten 
and producing at least 0.10 in. of runoff. The data for these storms are s 
in Table 2. 

The unit graph peak and time from beginning of rainfall excess to peal 
of runoff, tp, were plotted against rainfall intensity as shown on Fig. 5. 
plotted data indicate two separate groupings, depending on whether the ; 
of high intensity came at the beginning or late in the storm period. The 
on this figure (fitted by inspection) show that both the peak of the in 
and the time to peak are dependent on rainfall intensity and storm patter! 

To determine if similar relations existed on the larger areas, a ser 
storms was selected, for the 290-acre watershed, which had the a 
characteristics: total storm runoff not less than 0.10 in.; reasonably un 
rainfall intensity for the period of excess; and a duration of rainfall exces 
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than 10 min, or more than the time to peak. The group of storms select- 
ad a duration of rainfall excess from 11 to 30 min with the exception of 
}44-min period on April 22, 1944. 
rouping of the data according to the duration of the rainfall producing the 
graph, and plotting similar to Fig. 5, indicated that duration has little 
ing on the results so long as it is less than the time to peak. The plotted 
S again showed two separate groupings depending on whether the period 
sh intensity occurred early or late in the storm period. A comparison of 
lottings for these two watersheds showed that, for the larger watershed, 
all intensity has less effect on the peak of the unit graph and the time to 
One might logically expect that the runoff hydrograph peak resulting 
rainfall of a given duration would be more closely related to the net rain- 
ntensity, that is, rainfall intensity minus rate of infiltration. Attempts to 
slate this net rainfall intensity with the peak of the unit graph and time to 
gave no improvement. This was probably due to the difficulty of making 
mable estimates of the infiltration rates for such short periods. 
‘om the graphof Fig. 5,a seriesof storms was selected to represent vari- 
ainfall intensities. Unit hydrographs for each of these storms are shown 
zy. 6. The actual total runoff, for the storms selected, varied from 0.10 to 
in., with the exception of the storm of May 27, 1938, which had a total 
fof 0.66 in. The original unit hydrograph definition by L. K. Sherman, 3 
ied for use on very large drainage areas states: 


or the same drainage area, however, there is a definite total flood 
riod corresponding to a given rainfall and all one-day rainfalls, re- 
rdiess of intensity, will give the same length of base of the hydro- 
aph.” 


data presented here, which applies to very small areas, the base of the 
raph increases as the rainfall intensity decreases. Perhaps this is be- 
, for rains of less than 15 min., at low intensities, the soil must be at a 
f moisture content to produce a minimum of 0.10 in. of runoff than for 
of higher intensities. Only a small percentage of the runoff in these 
S was due to base flowor interflow since therate 3 hr after the peak was 
bout 1% of the peak and, therefore, no attempt was made to separate the 
zraphs into their component parts. 

veloping Synthetic Unit Hydrographs.—With a limited amount of data 
ble and with considerable scatter of points, as shown in Fig. 5, one may 
ntly find itimpossible to pick a storm representative of a desired inten- 
r storm pattern. 

tudy was made of the unit graphs to show times before and after the peak 
various percentages of the peak rate occurred. The data for watershed 
-e shownin Table 3. Similar data from the other two Edwardsville water - 
were included in this analysis. The study showed two distinct groupings 
he times were compared with rainfall intensity. However, all points fell 
group when plotted against the unit hydrograph peak. A plotting of the 
of occurrence of rates equal to 60% of the peak is shown in Fig. 7, in 
- indicates minutes before peak and + indicates minutes after peak. The 
3 show that, when times are plotted against the unit graph peak, one 
san adequately define the shape of the recession for all three watersheds. 
treamflow from Rainfall by Unit-Graph Method,” by L. K. Sherman, Engineering 
tecord, Vol. 108, 1932, pp. 501-505, 5 
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‘imes before the peak of the unit graph, being affected by rainfall intensity 
storm pattern, are not as consistent as times after peak, which are the re- 
of depletion of detention volumes after cessation of rainfall. The average 
) of times before to times after the peak, for various percentages of the 
, was determined for each of the watersheds. The results, shown in Fig. 8, 
cate considerable variation between watersheds. However, since less than 
of the runoff occurs before the peak, the use of an average curve should 
introduce serious errors. 

he data and analysis thus far is confined to those areas on which runoff 
surements are available. The results will have much wider application by 
ting them to ungaged areas. Points representing the unit graph peak and 
to peak for rainfall intensities of 1, 2, 3, 4, and 5 in. per hr, taken from 
surve of Fig. 5, and similar plottings on other watersheds at Edwardsville 
a 189-acre area near Hamilton, Ohio, are shown on Fig. 9. These curves 
y tothose cases where the period of high intensity occurs late in the storm. 
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20 


0 
0.20 0.30 0.40 0.50 0.60 0.70 
Average ratio of time before to time after peak 


FIG. 8.—RATIO OF TIMES TO PEAK 


‘urves show the effect of size of areaon the unit hydrograph peak and time 
ak. A similar set of curves can be prepared for storms having the in- 
: period at the beginning. The curves for times after the peak, at which 
us percentages of the peak rate was reached, were taken from plottings 
lar to those in Fig. 7. The curves showing times before the peak were de- 
ined by applying the average ratios from Fig. 8 to the times after the peak. 
series of synthetic unit hydrographs for a particular size area and differ - 
ainfall intensities can be obtained from Fig. 9. 

) illustrate this, assume a unit hydrograph is required for a rainfall in- 
ty of 4 in. per hr and a drainage area of 200 acres. To develop the syn- 
> unit hydrograph, enter Fig. 9 at 200 acres and go vertically to the rain- 
ntensity of 4 in. per hr, then horizontally from these points to find a peak 
of 1.85 in. per hr and time to peak of 22 min after the beginning of rain- 
xcess. Proceed horizontally along the line for a unit graph peak of 1.85 
er hr to read times before and after the peak at the percentages shown. 
:. 
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Estimated runoff, 1.81 in. 
Actual runoff, 1.79 in. 
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Estimated runoff, 0.51 in. 
J Actual runoff, 0.54 in. 
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Thus a rate of 60% of the peak would occur 5 min before and 12 min after | 
peak. 


COMPARISON OF COMPUTED AND OBSERVED HYDROGRAPHS 


The final step was to apply the time-retention curves and the unit hyd 
graphs previously developed to the actual mass rainfall diagram of a num} 


TABLE 4.—COMPARISON OF OBSERVED AND COMPUTED PEAK DISCHARGES 


March 


May 2.64 

July 1.00 

1939 July .63 

August .88 

1940 August .84 

1941 September .60 

1942 July 1.86 

December 43 

1,24 

1944 April 1,34 

August .30 

1945 June 1.40 

1946 August 2.87 

1,41 

1947 April 1.60 

1948 June 37 

July 50 

Bret 

1949 May 1,56 

June th 04 

1.83 

August 92 

1950 January 54 

.65 

1951 June 97 

sully ‘92 

<n 

87 

1952 March 1.45 
2.16 =a 

duly 1.89 

1953 April "73 

June 

1954 July te, 

September ‘430 


October 


of storms on both the 27.2-acre and 290-acre areas to see how the com 
hydrographs agreed with the actual record. The method used is illustrate 
Fig. 10(a) for the storm of June 16, 1945 on the 27.2-acre area (W-1). 
mass rainfall diagram shows considerable variation in rainfall intensity du 
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excessive part of the storm. The antecedent precipitation index was com- 
ed and used to select a time-retention curve, to be plotted directly below 
mass rainfall diagram. The vertical distance between the mass rainfall 
| the retention curve represents the potential runoff to be expected from 
nfall prior to that time. The difference in these ordinates at selected points 
es the increment of runoff to be applied to a unit graph selected from Fig. 6. 
The first trial was made using a series of 12 min intensities even though 
breaks in the rainfall diagram indicated considerable variation in intensity 
ing certain of these periods. Investigations showed that much better agree- 
nt with the actual hydrograph could be obtained by using variable time in- 
vals such that the intensity during any interval was more nearly uniform. 
> time intervals selected, as indicated on the mass rainfall diagram, varied 
nultiples of two and ranged from 4 to 12 min. For each of these intervals, 
unit graph was selected for which the intensity most nearly matched that 
he mass rainfall diagram. For example, the intensity from 2:11 to 2:33 was 
0 in. per hr and the unit graphused for this period was the one based on the 
rm of April 17, 1941. 

Using the same procedure but applying a series of synthetic unit hydro- 
phs obtained from Fig. 9 to thestorms of July 8-9, 1942 and June 9-10; 1952 
vatershed W-4 gave the comparisons shown in Fig. 10(b) and 10(c). Note that 
retention curves for each of the periods of excess is based on the API at 
beginning of that period. 

The final test of the principles involved was the computation of peak rates 
one or more of the outstanding storms each year on watershed W-1. Table 4 
ws that the difference between observed and computed peak rates is less 
1 10% for one-half of the storms. The greatest disagreement between the 
erved and computed hydrographs was for those storms beginning with a 
rt period of very high intensity. A slight difference in the position of the 
eretention curve makes considerable change in the estimatedrunoff for this 
ial period. The high rainfall intensity requires use of the unit graph with 
1ort, sharp peak, similar to the one derived for the May 27, 1938 storm on 
ershed W-1, which greatly magnifies any small discrepancy in the estimated 
off. 


CONCLUSIONS 


. There appears to be a goodrelationship, for this soil type, between ante- 
2nt precipitation and the amount of precipitation retained in a given time. 
'. Runoff from the most important storms can be estimated with reason- 
- accuracy from the mass rainfall diagram and time-retention curves. 

. These investigations show that on small watersheds, one unit hydrograph 
not adequately define the shape of the hydrograph derived from a storm of 
duration. 

For the small watersheds involved in these studies, there appears to be 
od relationship between rainfall intensity and the peak rate of the unit hy- 
sraph and the time from beginning of excess rainfall to peak rate of runoff. 
. A method has been prdsented for constructing synthetic unit hydrographs 
lrainage areas from 20 acres up to about 500 acres and for different rain- 
intensities. 

. Computed hydrographs show close agreement with the observed record 

e unit hydrographs used are based on storms having similar rainfall in- 
ities and distribution patterns. 


> 
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7. There is need for further research along these lines, particularly to a 
if similar relations exist for other soil types. . = 
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TWO METHODS TO COMPUTE WATER SURFACE PROFILES 


Discussion by S. V. Chitale and K. Arunachalam, F. Paderi 


Closure by Joe M. Lara and Kenneth B. Schroeder 


V. CHITALE! and K. ARUNACHALAM.2—The authors are tobe congrat- 
d for the lucid exposition they have given of the step-by-step method 
ted in the Bureau of Reclamation for backwater computations when both 
hanneland overbank discharges are required to be considered. Two meth- 
have been illustrated, Method B having shown to be more accurate of the 

Two assumptions have been made in the development of method; one is 
hg is the same for channel and overbank flow in a given cross section, and 
ther that velocity head is alsothe samefor both channel and overbank flow 
cross section. In some cases, these assumptions may not appear to be 

and their justification would depend on the magnitude of error involved 
ese assumptions. Application of another step-by-step method without re- 
se to the above assumptions has been illustrated below. It is contended 
this method is more general than Method Band its adoption in preference 
ethod B may be found desirable in specific cases. For comparison, an 
ple is solved by both the methods; that is, the general method as well as 
od B (Table 3). Since the data of Red Fox River given by the authors are 
sufficient to illustrate the more general method advocated below, another 
-—Brahmani River in the State of Assam, India—has been chosen for work- 
ut comparative results. 
neral Method 
ssumption.—The only assumption made is that water level across a cross 
on is same over the channel and overbank portions. This assumption is 
non to this as well as Method B. 
yuations.—Besides those given by the authors, the following additional 
ions are required to be adopted: 
sing the subscripts c and s to distinguish properties of channel and over- 
portion, at any section, 


.pril, 1959, by Joe M. Lara and Kenneth B. Schroeder. 
thf, Research Officer, Central Water and Power Research Sta., Poona, India, 
‘esearch Asst., Central Water and Power Research Sta., Poona, India. 
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h Veg - h vey = Ah ve see wales te we ee 
h Vgy ~ h Vg, = Ah Vg a toate tehedu de hwlaaPea ented i ome 


he, = AH + Ah vc - eddy losS ......-+e2eee 
and 


tl 


he, = AH + Ah vg - eddy loss ......-+ eee 


Ss 


Design Data.— Backwater computations were carried out for the Brahm 
River in Assam State, India: 


1. A plan map in Fig. 1 of the river shows the location and designatior 
the different cross sections. 


Petes. 7 oA al Pie Sa ee 


OVERBANK PORTION 


EXTENT OF 
---~--OVERBANK ---{------ 


i 


FIG, 1 


2. The water surface elevation at Cross Section 0-0 is 682.00. 

3. Discharge in the channel portion of cross section 0-0 at water surf 
El. 682 is 835,000 cfs. 

4. Discharge in the overbank portion of cross section 0-0 at water curl 
El. oe is 40,000 cfs. 

5. = 0.03 for the main channel and “n” = 0.05 for the overbank secti 

6. Values of Kd, A, R, etc., for different water levels are calculate 
given in Table 1. 

7. The relationships between water surface elevations and k dg for aitter 
cross sections are shown graphically in Fig. 2. 


Procedure. ; 7 

Step 1.—To start with, discharges occurring in the overbank and cha 
portions are known at the end section, that is, at cross section 0-0. In the f 
trial for estimating the discharges passing through the two portions at cI 
section 8-8, the velocity head is ignored, that is, hg = AH. Since AH isys: 
for both the portions, hg overbank = h¢ channel. A trial AH is assumed, Co 
Table 2. 
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p 2.—The values of Kd for the two sections 0-0 and 8-8 are averaged 
e discharges Q, and Qs are computed by Eqs. 2 and 3 and entered in 
10 and 17. 

9 3.—If the sum of Q, and Qg is not equal tothe actual total discharge of 
0 cfs, the actual discharge is divided in the ratio of Qc. and Qg obtained 
) 2 and adopted (Cols. 19 and 20) for finding the velocity head in the next 


TABLE 1.—BRAHMANI RIVER 


7 


Reach Bankful 
Elev, A P R Kd length hv Stage 


(a) Properties of Channel Portion 


39,665, 756 2.8174 
42,509,164 2.5880 
45,479, 809 2.3837 


28,546,541 
31,601,720 
33,944,904 


39,219,254 
41,901,816 
44,604,793 


46,177,066 
49,566,117 
52,990,611 


815,000 
1,910,000 
3,330,000 


1,470,000 
2,870,000 
4,630,000 


1,950,000 
3,430,000 
5,420,000 


1,340,000 
2,620,000 
4,240,000 


4.—Again a trial AH is assumed. This is added to the water surface 
n at cross section 0-0 to get water surface elevation at cross section 
¥ : 


5.—The velocity heads h ve and h vg using the above water surface ele- 


— by using the discharges obtained in step 3 (Cols. 4 and 11). 
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TABLE 2.—BACKWATER COMPUTATIONS FOR BRAHMANI RIVER Uj 


Cross|Trial |Water 
sec- | AH 


683,50] 4.50 
683,60} 4.50 


1,655}-0.395] 0.04 


Water 
surface 
elevation 


Section 


Trial II 684,00 


Trial III 684,15 


13-13 686.00 


Trail II 687.00 
Trial III 687.30 
Trial IV 687,34 
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%AL METHOD SUGGESTED BY S, V, CHITALE AND K, ARUNACHALAM 


Qs Q 2 


x 105 cfs 


Velocity head 
neglected 
2.67 

af 


Velocity head 
neglected 

1.30 ue 

1.60 1.09 

1,18 

1 


Velocity head 
neglected 
0.90 0.57 


0.58 
2 5 & 
eo V“Q x 10 hv |hvy - hvo 


0.57 
ae : 


Eddy} Mean 


Total 
loss 
1,320.00 
1,34 


16 


: 2,340.00 
0.67 5,60 | 4.17 0.90 [2.99 |3.89 | 2.09 
7.60 1,790.00 
1.15 17.70 | 3.22] 0.85 |0.43|2.59 |3.02 |2.17 
7.57 1,790.00 


17,70 0.43 |2.57 |3.00 


7.95 1,380.00 
0.80 5.20 2,98 
7.78 1,230.00 
0.97 7.10 3,21 
7.80 1,230.00 — 
0.95 7.10 3,19 
7.80 1,230.00 
0.95 7.10 3,19 
—~8 
“wy 
8 z on 
0 1,54 0.66 0.07 0.713 ]0.783 |1. 
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Step 6.—The terms Ah v, and Ah vg are calculated by Eqs. 5 and 6 
entered in Cols. 5 and 12. 

Step 7.—Eddy loss is assumed at 50% of Ah v for retarded flow and 1( 
Ah v for accelerated flow; Cols. 6 and 13. 

Step 8.—The values of hg, and hf, are then computed by Eqs. 7 and § 
entered in Cols. 7 and 14. 

Step 9.—The value of K d for the two sections 0-0 and 8-8 are averagec 
discharges Q, and Qs are computed. 

Step 10.—The discharges obtained in Step 9 still do not satisfy Eq. 1. 
get an idea of hg to be used in the next step, the following equation is tu 


woe geet 

Q4=Kd Lect 
Q . h 
actual Zz f ; 

rE LRT ee FETS approximately ..... a 

computed {computed 

Solving for hg: 
h -( Qactual 

Gocmpaieds f computed ot sina, aoe se Oe 


This assumes that Kd remains constant for small change in water sur 
elevations. 


ELEV. IN FT. 
ELEV. IN FT. 


30 
Kd IN MILLIONS 


40 
Kd IN MILLIONS OVERBANK PORTION 
CHANNEL PORTION 


7 
FIG, 2 


As the flow in the channel usually predominates, hg used in Eq. 11 
for channel portion. The terms Qactual 29d Qeomputed are the total 
charges. The value of hg obtained by Eq. 11 is Bakared: in Col. 21 for com 
son with the assumed hg, Col. 7 so as to give an idea of AH to be used i 
next step. oats a 

Step 11.—As we proceed, the discharges Qs and Q¢ approach the c 
value. This procedure is repeated until Eq. 1 is satisfied. The same 
dure is used to compute the upstream water surface elevations for cros: 
tions 13-13 and 16-16. 
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. PADERI.3—The writer has had occasion to compute water surface pro- 
for flow in open channels. He particularly appreciated the clear and ac- 
te exposition on the two methods of finite differences derived and is pleased 
them for several reasons, among them being: 


The development of Methods A and B is based ona profound and exten- 
analysis of the various physical phenomena involved and the relative ef- 
of each in the theoretical analysis. 

The authors point out that practical judgment, in addition to theoretical 
sis, also plays an important part in the evaluation of these phenomena. 

The writer noted in this paper several ideas he developed when dealing 
the same subject.4,5 


writer has nothing to add to the first two above statements, but some ex- 
tion of the third follows: 


A fundamental element of computing by finite differences as presented 
e authors is the hypothesis of linear variation of the function S¢ on each 
vision, specifically its equivalent 

Q2 
Ao XR 


e X, defined by Manning has a value in metric units: 


|= 


n American units: 


nie DARSR a ieingier Say stein ine N13) 


in metric units: 


1 American units: 2 
Srey pace Elsen ame yi ioeeh Saar Seis 9 6) 
1.486 A R® 


inear variation hypothesis, on each subdivision, of the ratio 


Q? 


A2x2R 


3st. Prof., Univ. of Pisa, Italy. 

sul moto permanente delle correnti a superficie libera,” by F. Paderi, Richerche 
sria, Roma, marzo-aprile, 1939, ’ 

Br oreit di rigurgito di correnti lente,” Stab. Tip. gia Chiari, Succ. C. Mori., 

e, 1946, e Bollettino Tecnico Ingegneridi Firenze, aprile-maggio, 1948, ~ 


e 


Nii 
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(not of each element such as A or,R, as used in previous studies) was alr 
stated, to the best of the writer’s knowledge, for the first time in the writ 
report of July, 1946, and presented at the Commission for the Professor 
(Chair of Hydraulics: D. M. June 28, 1946, and August 6, 1946). 

The above-mentioned hypothesis was used in the same report (applying 
method of finite differences) for the computation of water surface profiles 
subcritical flowin cylindrical channels. An example was shown in Section 
the report applicable to specific conditions. 

The writer appreciates the authors’ confirmation of the validity and p: 
ticality of the forementioned hypothesis, implicit in the evaluation of hg, 
each subdivision, with 


he = (L + 2) (Sty - Sto) J 0-614 ben ais ae 


2. A second element as important to practical hydraulic analysis as the 
discussed previously concerns the proper introduction in the computatior 
a Simplified and practical concept of eddy loss incident to the conversic 
kinetic to static head and vice versa in open channel flow. This is usually ¢ 
sidered as some portion of the change in velocity heads. 

In previous work4,5 (dealing mostly with subcritical flow) the writer 
veloped this same concept which was discussed in general terms. Base 
previous studies of this concept, the writer introduced in his conclus 
another coefficient, p (analogous to the authors’ coefficient C) combined: 
the coefficient X of uniform flow to evaluate the friction head and eddy lo: 
holding to the true phenomena. Regarding the choice of the p value whe 
constant or variable according to the type and condition of flow, the wr 
pointed out evidence for the necessity of conducting further research in part 
lar the carrying out of proper experimental tests. These tests should be gex 
to attain a better evaluation, both quantitatively and qualitatively, of the p 
efficient. 

The writer also explored the results of particular values of p, use 
practical applications ranging from zero to one and reciprocal acne 
the water surface profiles. For specific conditions, interesting evaluation 
the coefficient p or C, also designated as -p or k, were obtained in later s 
ies®,7 on this important subject. | 

Regarding these studies and those of the authors, it is important to note 
minute numerical departures are not bound to a particular single set of co 
tions which can be said of the general agreement of the fundamental phys 
concepts involved. Under certain circumstances the value of p, C, etc., ¢ 
also possibly be accounted for in the value of X itself which has n as the vi 
able factor. : 

It is well to note also that with the introduction in Bernoulli’s energy e 
tion of the eddy losses affected by the coefficient p, C, -p, or k, there rem 
in the said equation a single multiplier 1 - p (or 1 - C, 1+ p, etc.) of the 
ference invelocity heads. Bernoulli’s equation as quoted inpaper can be ¥ 
ten as follows considering only the eddy losses mentioned above: 


23.4 do 5 21.+ dy + bg +(1 -.C) (bys = hye) a 


6 “A duzzasztdsi gérbek szdmitasara ajanlott médszerek hidromechanikai 8sszel 
onlitdsa,” by Gy. Kovacs, Vi2tigyi Kézlemények, N, I, Budapest, 1952, pp. 84-92, 
“Etude et tracé des écoulements permanents en canaux et rivierés,” by R. Si 
Dunod, Paris, 1954, p. 13. s 
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‘ad of the factor (1 - C), a value of (1 + C)for this coefficient can be used, 
-an algebraic change in sign is noted. Based on the 100% value reported 
e authors for retarded flow conditions, the factor is resolved as follows: 


Na re a eS i I gna 6 ine mnie (18) 


authors are commended on their calling attention to the numerical value 
which can change to fit the conditions observed and its assigned value is 
‘(o the judgment of the hydraulic engineer. 


OE M. LARA8 and KENNETH B. SCHROEDER.?—The writers appreciate 
nteresting discussions submitted on the methods of computing water sur- 
profiles. 

lessrs. Golding, Anand, and Silvester introduced the idea of expressing 
eyance as K equal to 1.486 A R2/3. Adjustment of the Kg-values used by 
riters canbe readily facilitated in cases involving variations in “n.” This 
articularly adaptable to conditions where the water surface profiles have 
observed for a series of discharges which results in a definite “n”-stage 
tionship. The writers found this idea very helpful in subsequent studies. 
ir. Golding’s additional suggestion of replacing Col. 5 of Table 2 with two 
mns is also useful when applying the preceding analysis. He suggests spe- 
ally tabulating the values of K = 1.486 A R2/3 in one columnand Kg = K/n 
e other. Mr. Golding’s graph plotting the conveyance “k” vs. depth indeed 
d be helpful in the case of an infinitely wide flood plain. The graphs pre- 
d by Mr. Silvester in nomographical form, wherein the solution of the 
tion, Q = (K/n) (S1/2) is included for various prismatic channels also 
d be helpful in certain cases. 

he refinement of using different flow travel paths was questioned by 
srs. Golding and Peterson of Method B. This method has not been tried in 
cient cases by the writers to arrive at a more conclusive appraisal. In 
ew cases tried, it was also found by the writers that computed elevations 
not significantly affected. The method merits though from the standpoint 
nsidering the existent hydraulic conditions realistically, yet conforming 
andard theoretical analysis. 

he writers acknowledge with interest the work carried out by Mr. Golding 
e study wherein the “n” values were increased based on observed data. 
oncluded that when the deflection angle exceeded 45 degrees, the “n” value 
fiably could be increased. Although Mr. Anand referred to the formula, 
loss = coeff. x 1/2 (hvy + hy9)s he did not offer any suggestions or guides 


signing a numerical value for the coefficient. This, as with other hydraulic 
mptions, undoubtedly is left to the judgment of the engineer. Z 
essrs. Anand and Peterson offer further discussion of the eddy loss analy- 
From actual hydraulic measurements, Mr. Anand verified his selection 
dy loss coefficients which probably supports his recommending that eddy 
3s be neglected for converging reaches but a 50% correction used for a 
‘ging reach. Mr. Peterson, on the other hand, suggests a 10% correction 
onverging reaches and from zero to 100% for diverging reaches. 

r. Peterson’s presentation on developing the coefficient aj, and his analy- 
f the continuity equation resulting in the formula, Q = Kag Sgl/ 2, are of 
2mic interest. His example of computing the water surface profiles is 
iydr, Engr., Bur. of Reclamation, Denver, Colo. 

iydr, Engr., Bur. of Reclamation, Denver, Colo. 
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another approach incorporating the application of this a correction and us 
the geometric mean of the Kgs. The results, however, are quite different fr 
those of the writers. About a 50% error in discharges is indicated. 

It is to be remembered that in the final analysis of any method, regardl 
of the mathematics involved in satisfying Bernoulli’s energy equation, the co 
puted results carry only the same degree of accuracy as that of the obser 
field data. The moment eddy losses and corrections for velocity head chan; 
are introduced in any of these computations, the element of judgment by 
individual engineer plays a significant role. The writers believe, therefo 
that greater emphasis be placed upon gathering field data which are jud; 
representative of the hydraulic conditions prevalent inthe reach under inves 
gation. Whatever mathematical analysis is applied to these data is at the d 
cretion of the engineer. Undoubtedly many of us have our “pet” methods. 
The following comments apply to the discussion submitted by Messrs. Chit 
and Arunachalam: 


It was not the writers’ intention to show Method B was more accurate t 
Method A. The primary reason for presenting the two methods was to sl 
computational procedures covering situations where the flow lengths betw 
sections were different in the main channel and overbank area. However, 
mentioned previously, the results, applying Method A and using only the f. 
length of the main channel, do not vary significantly from those using Met! 
B as applied to various problems thus far. 

There are severalcomments regarding the supporting data that require « 
planation or clarification as listed below: 


a. There appears to be an error in the reach lengths quoted in Table 
They are listed as being identical for the main channel and overbank porti 
Reach lengths shown in Fig. 1 do not seem representative of the flow path tal 
by the water in the overbank areas. It seems more reasonable that the ff 
would follow more of sinuous path somewhat parallel to the main channel 
shorter. | 

b. It would be interesting to know how the hydraulic properties of e; 
cross section (Figs. 1 and 2, Table 1) for the overbank portion were det 
mined. The important question is: What happens when the bankfull stage 
reached and then surpassed? Normally, the cross section is laid out in 
field in a straight line extending to the limit of the overbank portion. Ift 
were done in the example, Sections 13-13 and 16-16 would intersect in 
overbank area imposing some complexity in the computation of the hydr 
characteristics. = 

c. A lesson in significant digits is to be gained inthe tabulation of the ve 
ous hydraulic characteristics listed in Table 1. In general, these can be | 
to three digits yet be compatible and representative of the gathered field 
A look at the tabulated data shows that they range from one place in R of 
overbank portion to eight places in the K dvalues for the main channel. T 
iS some inconsistency. : 

d. Velocity head (hv) values appear in Table 1 for all elevations of the 
channel of the four cross sections. No mention is made as to how or why 
values were tabulated. A check shows they were all computed on a basis 
discharge equal to about 855,000 cfs. 4 

e. In Table 1 there is also listed values for a term R.L. at each sect 
It is presumed this means ‘reference level.” However, to list them as b: 
full stage is another unexplainable item. 
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1 exception to some of the above comments, Messrs. Chitale’s and Aruna- 
lam’s method is judged meritorious and appears to be a worthwhile proce- 
>. The results of the discussers’ and writers’ procedures are plotted in 
_3 to show comparative differences. 

Mir. Paderi presents an interesting discussion of the results he attained in 
Studies involving finite differences. The writers have not had the oppor- 
ty to review any of the literature quoted by Mr. Paderi. It is noted, how- 
r, that most of his experiments and studies dealt with subcritical flow which 


KeE ioXtetoTHeINcn 359-5 
KEUFFEL@ESSERCO,  sective. ea 


4 5 6 
DISTANCE IN 1000 FEET 


FIG, 3 


writers believe he meant flows in the “shooting” range or at supercritical 
cities. Mr. Paderi’s comments relative to the coefficient assigned for 
ing eddy loss corrections are covered by previous discussors. The writ- 
appreciate Mr. Paderi’s remarks on the opinion he shares roparding the 
that there is still a lot of responsibility riding on the engineer’s shoulders 
e selection of these coefficients andit might be added that this applies also 
cement of “n” values for use in Manning’s formula. 
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FLOODS OF THE FLORIDA EVERGLADES® 


Closure by Edwin W. Eden, Jr. 


EDWIN W. EDEN, JR.1~Any discussion of the Everglades is incomplete 
ithout information on the Seminoles—its natural inhabitants—and their role 
the development of this unique area. There can be no doubt that the Semi- 
les have learned to live in the extremely wet, hot, and humid environment. 
ne flood problems arose when development by the “white man” required re- 
ef from the extremely frequent and long-duration flooding. Mr. Burnet’s 
scussion is greatly appreciated as it provides background on the Seminoles 
an extremely interesting phase—the development of the Everglades. 


hd 
z 
i 
- 


June, 1959, by Edwin W, Eden, Jr. : 
'Chf, Planning and Reports Branch, Engrg. Div., U.S. Army Engrg. Dist., Jackson- 
Corps of Engrs., Jacksonville, Fla, 
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the use of Eq. 8 for the solution of K. Since this is an analytical expressic 
the value of K can be solved to any desired accuracy. However, Eq. 8 is 
rather involved mathematical expression, a solution for K is not an easy matte 
It is believed that the graphical methods will provide a much easier and fast 
solution for those who are not mathematically inclined. For the sake of greate 
accuracy in using the graphical method, a probability paper of large sc: 
should be used. 

The study made by Mr. Sammons on the accuracy of the graphical metho 
is very interesting. When either the original variate or the log-transform 
variate is used, however, the result should be theoretically the same.2 A 
difference between the two approaches should resort to the personal error i 
volved in the graphical construction. 


2 Closure to “The Probability Law and Its Engi i ; 
; gineering Applications,” by V Cho 
Proceedings, ASCE, Vol, 82, 1956, pp. 6-7, Stele ee 


DETERMINATION OF HYDROLOGIC FREQUENCY FACTOR® 


Closure by Ven Te Chow 


VEN TE CHow, 1! M. ASCE.—The writer appreciates the constructive and 
‘resting discussions by Messrs. Weiss and Sammons. Both gentlemen sug- 
t the methods using the transformed variate. Mr. Weiss suggests an ana- 
cal approach, while Mr. Sammons suggests a graphical approach. 

It is believed that each method has its own merits and it should best apply 
1 case as required by the circumstances under consideration. The analytical 
thod requires only the direct solution of a mathematical equation. It can, 
refore, produce results to any desired accuracy, but the solution of the im- 
sit equation for the frequency factor is not an easy matter. The graphical 
thods replace the mathematical solution of an equation by the graphical con- 
uction. Their accuracy depends, of course, on the scale of the probability 
er used in the graphical construction. The use of an ordinary probability 
er should be satisfactory for practical purposes. Fora theoretical study 
ich requires high degree of accuracy, the analytical method is recommended. 
When y is normally distributed, its mean y is equal to its median My or 


ich is the value of y at a probability of “50%-of-time.” Since M is the 
dian of x, it follows 


My =InM....... catback aoe ete eee (2) 
om Eqs. i and 2, 

Y= INM .-cseesvvoes St Soa (3) 

M = exp.) oc odnewine a ere es Hs (4) 
stituting this expression for M in Eq. 6 in the original paper and solving 
f XK = exp (¥ + 02/2) ebn th ci eer ee Tae (5) 


Mr. Weiss uses this equation to derive an expression for a reduced vari- 
et. The latter is essentially a frequency factor in terms of the variable y. 
om Mr. Weiss’ Eqs. 3 and 4, the following may be written: 


: y/y = 1+C, I eS (6) 
i’ Oi oat ca wa ieere a eS (7) 
- Ky = t = oy/2 +(1/oy) In(1 + CvK)..... .(8) 


mps ring with Eq. 1 in the paper, it may be noted that Eq. 6 is an application 
ow’s general formula to the transformed variable y. Mr. Weiss suggests 


Q 


uly, 1959, by Ven Te Chow. ; 
rof, of Hydr. Engrg. in-charge, Univ. of Illinois, Urbana, Il. 
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NERALIZED DISTRIBUTION NETWORK HEAD LOSS CHARACTERISTICS 
Discussion by G. C. Anderson, J. V. Radziul, and P. Celenza 


§. C. ANDERSON, !6 M. ASCE.—The author has recognized that the method 
proportional loads is not applicable for determining the ability of a distri- 
ion system to supply fire-flow demands. After an arterial system, capable 
satisfactorily supplying the expected domestic and industrial consumption 
the design period, has been determined, it will be necessary to consider 
ividually all possible fire-flow demands on the system. This should be 
€ considering consumption at the maximum daily rate (the average rate of 
sumption on the maximum day) for the design period. The fire flow will 
e to be available from the arterial system at sufficient pressure so that it 
y be delivered through the secondary feeders and the minor distributers to 
desired fire location with a residual pressure of 20 psi. 


J. V. RADZIUL,17 and P. CELENZA.18_The author is to be congratulated 
his clear, concise documentation and presentation of a new method of cal- 
ation for the determination of overall system head losses through a wide 
ge of demand and equalizing storage rates utilizing a minimum of network 
lyses. His step-by-step solution of actual sample problems should be of 
istance not only to the hydraulic engineer in the design office but to the 
er works operator as well. 

it should be noted that some water distribution network designers may not 
sider that the functional relations of Eq. 1 represent an entirely new con- 
t,19 since various types of approximation techniques somewhat similar in 
ire have been used for years, in network analyses, to obtain crude esti- 
es of source to point system head losses under varying conditions of sys- 
. input without complete network calculations. However, this paper for the 
t time has proven that the functional relations of Eq. 1 can be satisfied only 
er proportional loading. Eq. 2 is a new and entirely original contribution 
he published literature. Application of the principles embodied in Eqs. 1 
2 willimprove water distribution network design and water system opera- 
a 

The writers believe that these procedures, in addition to augmenting design 
yputation techniques, provide a simple practical tool which can be readily 
lied by most water works operators or engineers in appraising the limita- 
iS of an existing water distribution system from either an engineering or 
rational viewpoint. The author’s contribution is indeed timely and valuable 


| January, 1960, by M. B. McPherson. 

6 Supervisory Engr., Natl. Bd, of Fire Underwriters, New York, N. Y. 

T Distribution Design Engr., Philadelphia Water Dept., Philadelphia, Pa. 

8 Asst. Distribution Design Engr., Philadelphia Water Dept., Phila., Pa. 

9 «Flow and Loss of Head in Distribution Systems,” by J. J. Doland, Journal, A. W. 


Vol. 35, No. 234, March, 1943. 
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when one considers that practically all of the approximately 18,000 water d 
tribution systems in these United States are experimenting “growing pains” 
are in the throes of substantial rehabilitation (Of the 18,000, about 300 : 
major systems.). 

With these procedures as an aid, the evaluation of an existing water dist 
bution system can be accomplished without considerable expenditure of ti 
or money for equipment. Good water works practice incorporates a yearly 
bi-yearly inspection of fire hydrants. Inspection procedure could easily 
amended to include a zero-flow pressure reading and the time of recordi 
This information, together with the elevation of the fire hydrant and norm: 
available pumping station and storage operating records, will establish su 
cient data for the use of the author’s equations. Computations can then 
made to determine head losses at a given point for any range of system dem 
and system input rates. 

Operators are confronted with the urgent need for pressure data, parti 
larly for the peak hour of maximum day consumption. This field data is di 
cult to obtain due to timing, coincidence of peak consumption loads and mz 
mum operation and construction work loads, personnel and equipment lim: 
tions, and the impossibility of saturating any sizeable area with press 
gages. The need is now obviated, since head losses for any given point car 
calculated directly for any system demand rate. 

Hydraulic grade contour maps for different levels of consumption car 
made for an existing system without recourse to a computer study or elabol 
hand calculations. The hydraulic grade contours superposed on a pipe sys 
diagram would provide a graphic picture of the system, indicating the stre1 
or weakness of the various areas. The study undoubtedly would confirm ¢ 
tain known low pressure locations and would probably reveal other conditi 
which are not readily apparent. Comparison of the slopes of hydraulic gradi¢ 
and pipe sizes would indicate those supply lines in which corresponding fl 
are much less than potential capacities. Overloaded supply lines would be 
dicated in the Same manner. Those locations for which the data is obviot 
inconsistent with calculated values obtained via Eqs.1 or 2 should bear furt 
investigation in the field. The writers’ experience indicates that inconsist 
cies can usually be traced to the influence of adverse local field abnormalit 
such as closed or partially closed valves, open boundary control valves 
tween pressure districts, or local non-proportional loads (severe oe 
of a sizeable industrial or commercial drawoff). 

In the author’s discussion on proportional load characteristics witl 
equalizing storage, he states with reference to Fig. 4 that “correlation is ¢ 
except for the “Southwest” station. Part of the scatter of field data, but 
probability of inexact proportionality of local loads must be considered at 1 
partially responsible.” The writers performed the network analyzer de: 
runs for the Belmont Gravity District shown in Fig. 4 and are in full agi 
ment with these statements. The “southwest” main is the principal fee 
serving a large industrial water consumer located approximately 5.9 air 
miles from the filtered water basin. During the period of field testing rey 
sented in Fig.4, pump operation at this industrial plant caused pressure f 
tuations of 3 to 10 psi; the resulting water hammer when the pump was sta 
or stopped produced a local short-time pressure fluctuation of upto 38 psi 
damped effects transmitted clear back to the filtered water basin. This s 
ation was remedied by having the pump draw directly from an existing sto! 
reservoir at the plant, filling the reservoir from the system. The writers 


: 
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ident that closer correlation would not be obtained with new field measure- 
ts at the “Southwest” point since “calculated” head losses to the plant (fed 
§-in. main shown at bottom of Fig. 3) are in good agreement with current 
ational data telemetered tothe Load Control Center, which includes a con- 
mus record of pressure at the plant. 

ig. 6A shows the proposed Belmont High Service District for the near fu- 
and incorporates proposed improvements consisting of a 20-in. reinforc- 
main and an elevated tank. The existing 1960 district is a closed pumping 
em without storage. Current data for this district covering a wide range 
ow data selected at random from the log sheets of the City of Philadelphia’s 
1 Control Center appear in Table 5. Elevation of the station discharge 


TABLE 5.—BELMONT HIGH SERVICE DISTRICT (1960) 


PR; Gage Dhj, in feet, from Pp to Point X 
Discharge} Qq = Qp hae Measured Calculated 
Pressure| in mgd | *) £0! “| (FieldData) [m= 1.85 | m = 2.00 
, in psi in psi 
| 3/14/60 53.8 6.4 64.5 17.0 17.3 16.8 
0000 | 
3/14/60 54.0 10.0 54.6 40.3 39.6 40.9 
} 1830 . 
3/14/60 SiG- 6) 38,6 63.0 26.28 26.2 26.2 
2330 
3/20/60 58.4 ; 8.5 64.0 28.8 29.3 29.6 
1630 
4/2/60 B35. 5h. 0.6.2 66.0 12.8 16.3 15.7 
| 0000 
| 4/2/60 48.4. |. 5.0 62.3 9.6 
‘| 0200 . 
| 4/2/60 54.0 9.2 57.2 34.3 
| 1300 
‘| 4/3/60 66.6 7.8 
| 0200 
| 4/3/60 } 63.4 22.7 
| 2330 


Used to determine K, for “calculated” head losses. For m = 1.85, K, = 0,559, For 
2,00, Kj = 0.409. 


sure gage at Pp is 258.3 ft. The pressure gage at Point X is located on 
-in. main, very near the site of the future elevated storage tank shown in 
6A. This 8-in. main comprises part of the distribution grid of a small 
area supported by the 12-in. main shown in Fig. 6A and an old 8-in. and 
d 12-in. feeder in the vicinity of the future 20-in. main. Elevation of the 
pressure gage at Point Xis 216.6 ft. Item No.3 in Table 5 was arbitrarily 
ed todetermine Kj. Since Qg = Qp at Pp, Eq. 1 applies. For Dhy = 26.2 
ith m taken as 1.85, Ky = 0.559; with m taken as 2.00, Ky = 0.409. The 
culated” head losses from Pp to Point X, when using either m = 1.85 or 
ompare quite favorably with the “measured.” With the exception of 
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Item No. 5 in Table 5 “calculated” values are within 1/2 lb of the “measur 
values or indeed well within the normal combined probably error of the gag 
The Belmont High Service District is predominately a residential system 
relatively minor industrial usage. The close agreement of the results 1 
actual values, for this system at least, lends further validity to the use of £ 
portional loading as a design assumption. Inconsistencies between an ear 
field survey and calculations made at that time by the author led to the loca 


TABLE 6.—TORRESDALE HIGH SERVICE AND FOX 
CHASE BOOSTER DISTRICT (1960) 


Dh, in feet, from ] 
to Point Y 


3/18/60 
1030 


3/18/60 
1630 


3/19/60 
1530 


3/19/60 
2200 


4/2/60 
1700 


4/3/60 
0300 


6/13/60 
1200 


6/13/60 
1800 


6/14/60 
0600 


6/14/60 
2000 


4 Used to determine ¢ and n for “calculated” head losses. For m= 1.85, = 0. 
and n = 2,60, For m= 2,00, = 0.421 and n = 2.55. 


of an 8-in. and two 6-in. boundary valves which had been inadvertently ope 
cauSing a leakage of about 15% of the HighService output to the Belmont Gr: 
District. 

Fig. 5 shows the future Torresdale High Service District with proposed 
ing improvements indicated by dashed lines. The existing (1960) distr 
also served by storage from Standpipes at A. It incorporates an mee, 
area, called the Fox Chase Booster District, as indicated on the lower ed 
Fig. 5. This present composite service district is supplied not only fror 
but also from either the Oak Lane Pumping Station or the Fox Chase Bo 
Station via the Fox Chase District (both pumping stations are located belo’ 
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stem shown in Fig. 5). Random data for the existing district, with the Fox 
ase Booster shut down, taken from the log sheets of the Load Control Center 
year in Table 6. Elevation of the station discharge gage at Pp is 0 ft. The 
sSure gage at Point Y is located on a 12-in. grid main in the vicinity of 
int B in Fig. 5. Elevation of the field pressure gage at Point Y is 154 ft. 
¢ values from Items 5 and 6 were substituted in Eq. 2. Simultaneous solu- 
n gave @ = 0.542 and n = 2.60 for Point Y, for m = 1.85; @ = 0.421 andn 
.55 for m = 2.00. Comparison of the “calculated” versus the “measured” 
id losses from P-y to Point Y when m = 1.85 is considered good and within 
septable limits of accuracy. Moreover, the Py discharge pressure was only 
ainable to the nearest point. Load Control Center personnel stated that the 
or in measurement from Py to Point Y could be as muchas 1.5 psi, depend- 
on operational circumstances. Torresdale High Service District is rela- 
ely residential in character with large wooded tracts and farm areas inter- 
srsed with localized concentrations of industry or institutional developments. 


TABLE 7.—HAZEN-WILLIAMS C - VALUES - EXISTING DISTRICTS 


Belmont High Service (Existing Arterial Piping in Fig. 6A) 


a oe 
20 


jiameter, inches 


-In Fox Chase District: 20 in., C = 124; 16 in., C = 48; 16 in., C= 44; 12 in.,C = 47. 
= 7 


re, in this district, utilizing equalizing storage, the limitations of propor - 
nal loading are apparently applicable as wellas the assumption that the pipes 
he average fit m = 1.85. <a 

he use of the author’s generalized network head loss characteristics in 
unction with field data has been demonstrated in this discussion. The equa- 
s under the limitations of assumed proportional loading and the application 
= 1.85 for “new” and “old” pipes regardless of their condition give “cal- 
ed” head loss values infairly good agreement to actual “measured” values 
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for the two Philadelphia districts cited. In the case of the Belmont High | 
vice District, use of m = 2.00 gives “calculated” values that compare as fa 
ably with the “measured” values as those obtained when using m = 1.8) 
= 1.9would provide the best fit between actualand calculated values. It sh 
be noted that the gage at Point X is at the extreme end of the Belmont . 
Service District and is located in a weak, local distribution grid which w 
be subject to more severe fluctuation than a point nearer the arterial pij 
For the Torresdale High Service District it was found that m = 1.85 prov 
the best fit as opposed to any other value for m. In the writers’ opinior 
use of m = 1.85 for network analyses is reasonable and consistent when no! 
design assumptions and overall accuracy of data and results is conside 
In Table 7 are listed the Hazen-Williams C-values for the Belmont High | 
vice and Torresdale High Service Districts. With such wide ranges 0 
nee it appears that m = 1.85 for “new” and “old” pipes is thoroughly ju 
fied.20 

For flow to storage the author states that Qp/ Qq > 1. Obviously this aly 
holds true only when a district with equalizing storage has only one pum 
station as illustrated in Fig. 5 and by Table 1. In districts with equali 
storage and multiple pumping stations Q,/ Qg > 1 only when the storage ir 
rate is greater than the combined output of the other pumping stations 
Items 4 and 6 in Table 6 and the No. 4 series of runs in Table 4). 

The foregoing suggests the use of the author’s concepts for operational 
poses. High and low “normal” operating pressure values could be determ 
for critical areas to obtain criteria for activating warning devices. Co: 
centers could establish valid operating standards which would result in 1 
efficient and economical operation. In complex systems using two or 1 
sources of input plus equalizing storage, more effective utilization of the 
tions to maintain adequate district pressures would be assured regardle: 
tank inflow or outflow rates. Conversely, improper station operation or ir 
ularities in the distribution system would also be indicated. 


+ 
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ulas,” Proceedings, ASCE, March, 1959, by M. B. McPherson, Proceedings, — . 
September, 1959, pp. 143- 185: 


THE FOURTH ROOT n-f DIAGRAM@ 
Discussion by Donald Van Sickle 


SONALD VAN SICKLE,23 A.M. ASCE.—There has been presented a “gen- 
lresistance diagram” for use with the Darcy-Weisbach formula (and with 
Manning formula as well), presumably as a substitute for the much- 
roduced, but apparently seldom-used, Moody diagram. The merits and de- 
rits of the “general resistance diagram,” as opposed to the common em- 
ical formulae, have been discussed frequently in the various technical jour- 
324,25, 26 
Since Mr. Blench is presenting the diagram as a practical aid and has pur- 
ely avoided detailed theoretical explanations, further discussion of the theo- 
ical points is beyond the scope of the paper. There are, however, a number 
eatures of the diagram which merit some discussion. 
First, the proposed diagram incorporates the Blasius equation to represent 
100th-turbulent” flow conditions, rather than the familiar von Karmdan- 
ndtl logarithmic relationship of the Moody diagram. Basis for this change 
he belief that the logarithmic formula lacks a “completely logical founda- 
1 and verification in the range of practical data.” The writer feels, how- 
r, that most of the published data appear to substantiate the von Karman- 
ndtl curve. The data presented by Peter Anton Lamont,27,28 with one ex- 
tion, give considerable support to the von Karman-Prandtl curve, and the 
eriments on newaluminum pipe,29 on a large, concrete-lined tunnel, 39 and 
coated steel penstocks,31 offer further evidence that the Blasius equation 
ot applicable “in the range of practical data,” except at Reynolds numbers 


1 January, 1960, by T. Blench. 

3 Hydr. Engr., Turner and Collie Cons, Engrs., Inc., Houston, Tex. 

24 «Flow in Rough Conduits,” by Henry M. Morris, Jr., Transactions, ASCE, Vol. 

, 1955, p. 373. 

5 “Design Methods for Flow in Rough Conduits,” Proceedings, ASCE, Vol, 85, July, 

p. 43. 

: . “Relationships Between Pipe Resistance Formulas,” by Walter L, Moore, Pro- 
ings, ASCE, Vol. 85, March, 1959, p. 25. 

IT “A Review of Pipe-Friction Data and Formulae, With a Proposed Set of Exponent- 

formulae Based on the Theory of Roughness,” by Peter Anton Lamont, Proceedings, 
1e Inst, of Civ. Engrs., Part Il, Vol. 3, April, 1954, p. 248. 

28 “Formulae for Pipe-Line Calculations—Technical Committee Report B,” Pro- 

= gs, of the Third Congress of the Internal. Water Supply Assn., London, 1955, 

9 “The Prediction of Flow Rates in Aluminum Pipe, Irrigation Tubing, and Fittings, 
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below 2 x 105. Above this limit, use of the Blasius equation will give f-va 
and friction losses lower than those which have been determined by exp 
ment. 

Second, Mr. Blench defines his “roughness height” in terms of a base w 
is different from the base for k in the Moody diagram. Most of the publi 
roughness-height values have been expressed in terms of the “equivalent 
grain diameter” of Nikuradse. It would seem as though any departure f 
this established practice would cause considerable confusion and can be jt 
fied only under extraordinary circumstances. The reason given, howeve 
to “give e-values that suggest something like the over-all height of irregul 
ties in unlined rock tunnels... in deference to... the preference. . 
values that suggest visual roughness heights.” 

This attempt to give some physical significance tothe e-value is, of cou 
commendable. Unfortunately, the practical value of such a move is quest 
able since, as the author notes, the ‘irregularities in unlined rock tunnels 
“the only practical roughnesses that are fairly common and are really ass 
able by eye.” For lined tunnels and all types of pipe conduits, then, the 
posed e-values will have no more physical significance than those on whicl 
Moody diagram is based. 

It should be noted that the author “would be content to treat e asa 
number to be attached to roughness types that would be recognized intuitivs 
Such a code number has, apparently, not prevented the Hazen-Williams 
Manning equations from gaining widespread acceptance. Perhaps Mr. Bl 
would be willing to reconsider this change in the definition of e to prevent 
fusion. 

Third, the proposed diagram eliminates the family of transition curves 
tween the “smooth turbulent” and “rough turbulent” flow regions which 
characteristic feature of the Moody diagram, and replaces them with a g 
of curves that represent the trends for different materials. The reasons ¢ 
for this change are that the family of curves represents extrapolation wel 
yond the range of Colebrook’s data, and that the transition curves are “to 
unlike ones that can be found for other materials.” A comparison of Fi 
and 2 reveals that Colebrook’s data covered a very considerable portion o 
transition region. Of the twenty transition curves shown in Fig. 1, onl 
upper two and the lower three seem to represent a significant extrapola 
Since there is a paucity of data in both of these areas, it is difficult to a 
the validity of these extrapolations, particularly in view of the abundan 
data corroborating the curves in the middle portion of the transition re 

Presumably, the curves mentioned for “other materials” are those 
Nikuradse’s experiments with uniform sand roughness, 32 and from the 
by G. H. Hickox, A. J. Peterka, and R. A. Elder33 with discussion, on | 
steel and concrete lined tunnels. As noted, these curves are unlike the M 
diagram transition curves. 

It is implied that the Colebrook transition curves should be used onl 
“commercial iron pipes,” and the other curves used for concrete pip 
large tunnels. Lamont’s data,27,28 however, show asbestos-cement — 
cement-lined pipe, and concrete pipe conforming to Colebrook’s trans 
curves. In addition, there does not appear to be any unanimity of opinion i 
literature that concrete surfaces are comparable to the Nikuradse uni 
eS ee a eee eee! as ee ee 


32 Engineering Hydraulics, by Hunter Rouse, John Wiley and Sons, Inc., New 
1950, p. 112. : 


33 “Friction Coefficients in a Large Tunnel,” by G. H, Hickox, A. J. Peterka, : 
A. Elder, Transactions, ASCE, Vol, 113, 1948, p. 1027. 
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id surfaces. To the writer’s knowledge, there are no reliable data available 
concrete surfaces which indicate that the Nikuradse transition curves are 
ly representative of these surfaces, and he must conclude that the Nikuradse 
"ves have been included simply to illustrate that some types of surface can- 
‘be represented by the Colebrook transitions. 

It is pointed out that the Hickox, Peterka, and Elder curves (and presumably 
Ontario tunnel curve at the bottom of Fig. 2) “show the kind of trend to be 
pected from the N curves, but to a different degree.” Hickox, Peterka, and 
ler in their discussion of these data conclude, however, that “although the 
alachia tests were carefully made, they do not conform to either the 
suradse or the Colebrook type of transition.” Cyril Frank Colebrook,34 in 
lition, in his discussion of the Ontario tunnel tests, comments that “despite 
appreciable experimental scatter, the test results are in very satisfactory 
“eement with” the von Karman-Prandtl smooth-pipe law. 

The apparently conflicting interpretations of the meaning of these curves 
d the writer to conclude that it would be less advisable to extrapolate from 
m to establish a curvefor a particular conduit than to extrapolate from the 
lebrook transition curves, for which a fairly large body of data is available. 
rtainly, before one could reasonably extrapolate from these new transition 
‘ves, it would need to be shown that they are truly representative of the 
lilable data on large conduits and tunnels. 

An examination of the twenty-nine curves shown in Fig. 19 of Hickox, 
terka, and Elder’s paper (not all of which are for large diameter conduits) 
icates that about half of the large conduit curves have the same trend as the 
lebrook transition curves, several are similar to the Nikuradse transition 
‘ye, and the six selected by Mr. Blench represent the most extreme depar- 
e from the Colebrook transition curve. It would seem to be difficult to de- 
‘ibe them as being “representative.” 

Colebrook has also presented additional data for lined, unlined, and partly- 
edrock tunnels.35 In the discussion of this paper36 are presented additional 
‘ves for large diameter conduits. The general trend of the curves, with the 
‘eption of the Apalachia results, follows the Colebrook transition curves, 
hough Colebrook cautions that “the experimental scatter is generally too 
sat... to determine the actual shape of the transition curve for concrete 
faces.” The data for the new Niagara tunnel30 fall in a pattern somewhat 
nilar to the curve for the Ontario tunnel shown in the author’s diagram. 
B. Bryce and R. A. Walker, however, make no attempt to relate these data 
either the Colebrook or Nikuradse transitions. It would appear appropriate 
Mr. Blench to include a number of these additional transition curves to 
ke it clear that the data for tunnels and large conduits are not allas differ- 
‘from the Colebrook transition as his initial selection would suggest. 
Finally, the proposed transition zone curves are said to be “widely spaced 
‘ause the practicaluser is expected to sketch transitions in terms of (a) the 
en ones, (b) allother available information relevant tothe specific problem, 
1 (c) personal judgment.” The wide spacing and completely different shapes 


eS ee 
34 «Turbulent Flow in Pipes with Particular Reference to the Transition Region Be- 
en the Smooth and Rough Pipe Laws,” by Cyril Frank Colebrook, Journal, Inst, of 
jee Vol. 11, No. 4, February, 1939, p. 133. 
35 “The Flow of Water in Unlined, Lined, and Partly Lined Rock Tunnels,” Proceed- 
, of the Inst. of Civ. Engrs., Vol. 11, September, 1958, p. 103. 

Discussion of “The Flow of Water in Unlined, Lined, and Partly Lined Rock Tun- 
. of the Inst, of Civ. Engrs., Vol. 12, April, 1959, p. 523, 
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of the curves, together with the usual lack of specific roughness informe 
available: in designing a new conduit suggest that “personal judgment” w 
be a critical factor in the use of Fig. 2. Since most designers have not 
Mr. Blench’s very considerable experience in conduit design, this import: 
of the judgment factor would seem to limit the usefulness of the diagram. 
use of additional representative transition curves would increase the diagré 
value to less-experienced designers. 

Of the objections to the use of the Moody diagram listed in Appendix 1. 
writer questions the validity of items i, ii, and iv, but must agree with it 
iii, v, and vi. In connection with item iii, however, Peter Ackers37,38 
published a set of design charts which considerably simplify the use ot 
Moody diagram. It should be pointed out, however, that the charts are b: 
on the Colebrook transition curve and would not be applicable for diffe 
shapes of transition curves. 


1 


37 “Resistance of Fluids Flowing in Channels and Pipes,” by Peter Ackers, H 
Research Paper No. 1, Her Majesty’s Stationery Office, London, 1958, 

38 “Charts for the Hydraulic Design of Channels and Pipes, ” Hydr, Researeh P: 
No, 2, Her Majesty’s Stationery Office, London, 1958, 


HYDROLOGIC STUDIES BY ELECTRONIC COMPUTERS IN TVA? 
Discussion by J. L. Kovner 


J. L. KOVNER.6—The use of automatic electronic computers in processing 
j analyzing hydrologic data is increasing ata rapid rate. Mr. Snyder has 
jwn one cost comparison between hand and machine methods which illus- 
tes why this development is taking place. In addition, it is pointed out that 
sk calculators have definite limitations in handling multiple regressions 
th ten independent variables. 

The author mentions other basic computations suchas mean daily discharge 
i sediment concentrations. These are not complicated calculations but the 
lume of data favors electronic computers. The USGS has a new water-level 
corder which punches head measurement at regular time intervals on paper 
ye and translates them to magnetic tape for immediate automatic processing. 
The general application of the machines in solving engineering and statistical 
oblems makes it possible for TVA to have its own installation. Those not 
fortunate can purchase machine time at normal rates from computer labo- 
tories. No one engaged in extensive hydrologic investigations need forego 
> computing improvements described by Mr. Snyder. 

The writer spent some time on the four watersheds prior to 1955, in the 
rly stages of the experiment. The design represented quite a departure 
ym the standard practice of using control watersheds. It was apparent, then, 
1t the four watersheds would have to be calibrated to account for varying 
steorological conditions since the 4 by 4 Latin square design would not ef- 
stively control these variables by the process of randomizing the treatments 
time. The writer was skeptical then that this could be done even though it 
peared that the small watersheds with uniform cover should behave almost 
runoff plots. Mr. Snyder admits that the old techniques did not work and 
may be said that the high speed computing machines are the modus operandi 
the current state of affairs. 

The author is certainly to be congratulated in developing this method of 
alysis. Although the final relationship for storm peaks is expressed explic- 
y in Eq. 5, one should not overlook the manner in which various hydrologic 
ocesses are taken into account, using a complicated bookkeeping system. 
is represents definite progress in hydrologic analysis and further applica- 
mn of the techniques is warranted. 

It appears, however, that the paper has some shortcomings in evaluating 
2 results. The experiment was not set up to test whether the cultural treat- 
ents would have any effect. There was enough evidence on hand to verify that 
west North Carolina row crops on steep slopes would produce much higher 


February, 1960, by Willard M. Snyder. 
6 Statistician, Rocky Mountain Forest and Range Experiment Sta., Forest Service, 
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storm peaks than pasture. The experiment could not demonstrate agai 
these facts even on the basis of low precision. The analysis must supply 
formation onthe difference between treatments within fairly narrow confide 
limits. 

In this respect, the results in Table 2 are not as good as indicated, nami 
that differences between peaks for various covers appear to be real; this is 
enough information. With standarderrors of 50% and higher—even for past 
cover—it will require great improvement in the regressions to obtain satisf 
tory confidence limits. For watershed No. 1, pasture, the 95% confide 
limits are approximately -0.07 to 0.93 cfs. The predicted peaks used dey 
considerably from the means in Table 1, and this accounts for the large st 
dard.errors. These predicted values are most interesting, however, si 
from a land-use point of view, it is the large storm peaks that are import 

In this connection, the use of the regressions for tests of significance 
confidence intervals as suggested is open to some question. The assump 
of normality of the dependent variable for a net of fixed values of the indepe 
ent variables would have to be verified. Storm peak data are usually sign 
cantly skewed from normal. It also appears from Table 1 that variances 
crease with size of peak. 

It is suggested, under the heading “Equation for Peak Discharge,” that 
watershed 5 special analysis will be needed to adjust for changes in cover 
is evident that the regressions for corn and wheat on any watershed will h 
larger variances than pasture because of changes in cover from bare soi 
planting time to partial plant cover at full establishment. 

In effect, separate regressions computed for the storms indifferent mol 
would have different parameters in Eq. 5, and a combined regression for 
storms may not be justified. The analysis could probably be improved by ci 
paring regressions on a Seasonal or monthly basis. 

Mr. Snyder uses an arbitrary and fixed time breakdown of the storm per 
There are four periods for storm peak hydrograph analysis and five per 
for storm runoff. These were evidently adopted on the basis of an ove 
analysis of runoff-producing storms. It would be interesting to know whet 
further subdivision of the storm duration into contributing periods would 
prove the relationships. The ultimate in electronic computing for this typ 
analysis appears to be a continuous flow of input data of the hydrologic 
meteorological variables. 


SCOUR AT BRIDGE CROSSINGS4 
Discussion by Joseph N. Bradley 


JOSEPH N. BRADLEY,13 M. ASCE.—The model results on scour at bridge 
iiments and piers have been treated in a logical and commendable manner. 
= consistency with which the model results plot could be misleading how- 
sr, by making the prediction of scour appear as a simple routine procedure. 
ce Mr. Laursen has said little with regard to the limitations of his results, 
ew remarks on this phase of the subject may be appropriate. 

Figs. 5 and 6 represent the results of model studies on abutment scour 
de under essentially ideal conditions, for example, with bed of granular ma- 
jal free to move, a rectangular flow cross section, and a uniform velocity 
tribution. One can compare the author’s results with those of a completely 
ependent set of experiments made with different size and gradation of bed 
terial, a constant depth of flow and a uniform velocity distribution, per- 
med at Colorado State University, 8 and find that the two are inclose agree- 
nt. This indicates that the model results are consistent and easy to dupli- 
e. 

Where streams meet similar specifications in nature as those of the models, 
re should be a reasonable correspondence between model and prototype. 
ere are streams in India and Pakistan which do approach these so-called 
al conditions; the stream beds present an unlimited depth of alluvial ma- 
ial, the river channels are extremely wide with more or less constant width- 
depth ratio, and the velocities are low due to an unusually flat gradient (the 
rage is 2 ft in 4 miles in East Pakistan). Under such conditions the velocity 
tribution cannot vary greatly across the stream. Field measurements of 
muir at bridge abutments, spurs, and guide banks for the rivers of India and 
<istan are on record!4 and these show surprisingly good correlation with 
‘model results.15 

Limited experience with scour on rivers in the United States show less fa- 
‘able comparisons. The reason is obvious, the gradients are steeper, the 
ss sections are irregular, the velocities are higher, and the velocity dis- 
butions are far from uniform. Under these conditions the greatest scour 
§ not necessarily occur at the abutments but is more likely to be found in 


4 February, 1960, by Emmett M, Laursen. 
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14 «The Behavior and Control of Rivers and Canals,” by Sir Claude Inglis, Research 
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5 “Field Verification of Model Tests on Flow Through Highway Bridges and Culverts,” 
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the portion of the channel where the depth of flow and velocity are greate 
Records of the United States Geological Survey of bridge sites in Mississig 
where the beds are generally of alluvial material, show this to be true. Th 
is also evidence from past floods in various sections of the United States t 
the settlement of piers in the deeper portion of the channel is more comn 
than abutment failures due to scour. Yet if Figs. 3 and 6 are consulted, 
order, it is found that prediction from the model gives scour depth up to th: 
times the pier width for the center of the channel, while scour up to six tin 
the depth of flowis supposedly possible at abutments. The latter can resuli 
fantastic figures, which are true in the case of the model with rectangu 
cross section, where allscour is concentrated at the abutments; but such p 
dictions are unreasonable when applied to irregular cross sections in the fie 

This discussion was not written to confuse the issue or discredit the moa 
results (which are valid for the conditions tested). Rather it is to point 
some of the remaining unknowns and (1) encourage investigators to mak 
concerted effort to take soundings of streams at constrictions both before | 
during floods for the purpose of better understanding the field problem 
making better utilization of the model results; and (2) towarn engineers to 
use the model information blindly but to treat each river crossing as an ir 
vidual problem, using the model results as a guide rather than a definite sc 
tion. Returning to item (1), further model studies will not produce the ans\ 
desired; field measurements are the only alternative. A further comment 
item (2) is that model results applied to abutment scour in the United Ste 
will certainly fall on the side of safety. The extra cost of unduly deep footi 
may, therefore, be sufficient to finance a comprehensive field study in a v 
short time. 


\P EFFICIENCY OF RESERVOIRS, DEBRIS BASINS, AND DEBRIS DAMS 
Discussion by Herman G. Heinemann 


iTERMAN G. HEINEMANN, !7 F. ASCE.—This paper is a very helpful sum- 
*"y of the better-known articles on the subject of trap efficiency. It also 
$ call attention to the need for more study, research, and analysis on this 
‘emely important matter. 

peaking generally, there are three sedimentation considerations of impor- 
e that face the planner or designer of a proposed reservoir: (1) How much 
iment will be delivered to the proposed reservoir (sediment yield)? (2) 
it percentage of the delivered sediment will be retained in the reservoir 
p efficiency)? (3) Where will the retained sediment be deposited (sediment 
ribution)? Help on any of these questions is a definite contribution to the 
wledge of reservoir sedimentation. 

t is assumed that the definition of trap efficiency given inthe “Introduction” 
not intended. Trap efficiency is generally expressed as the ratio between 
ment accumulation and sediment inflow. 

‘or clarity purposes, it should also be pointed out that the term “sediment 
juction” as it is used under the heading “Summary of Paper by Brune” and 
. in the referenced article by Glymph, has been largely replaced in the 
recent years by “sediment yield.” This also pertains to the heading of 
lin Table 1. This is made apparent by the use of the term “sediment 
d” in the later portions of this paper where reference is not made to arti- 
or papers of a number of years ago. These terms are still used somewhat 
rchangeably, but sediment production is now considered as referring to 
sion or sediment volume from small areas—plots—single cropped, etc., 
reas sediment yield refers to the sediment volume from more complex 
srsheds or subwatersheds—as defined by Glymph.18 

ince the trap efficiency for a structure is a changing value, high at first 
lower as the reservoir becomes filled with sediment, it is the writer’s 
ion that a date or a period of time should always accompany the trap- 
‘iency value given for any reservoir; it was expressed in this way for a 
ber of reservoirs included in the paper. One should not give a trap effi- 
ey value for “the life of a structure” without defining this or giving the 
od of years considered as the usefulor design life. Trap-efficiency values 
be more meaningful if expressed as follows: 

“The 50 yr. trap efficiency of x reservoir is estimated at Y%.” 


2) “The trap efficiency of x reservoir is Y% (1959).” 


| age 1960, by Charles M. Moore, Walter J. Wood and Graham W, Renfro. 

‘Hydr. Engr., Sedimentation Studies, Watershed Tech, Research Branch, Soil and 

Conservation Research Div., Agric. Research Service, Lincoln, Nebr. 

3 “Studies of sediment yields from watersheds,” by Louis M. Glymp, Jr., Internatl, 
ff Geodesy and Geophysics, Tenth Genl. Assembly, Rome, Italy, 1954, Vol. 1, pp. 

91, 1955, illus, 4 
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TOLKMITT’S BACKWATER AND DRAWDOWN CURVE TABLES? 
Discussion by Ven Te Chow 


VEN TE CHOw,!9 F. ASCE.—Tolkmitt developed his tables on the assump- 
1 of a parabolic channel having the width sufficiently large compared with 
depth y, so that the wetted perimeter P can be considered approximately 
al to the top width T. Accordingly, T varies with y0.5; the water area A 
ies with y1-5; and the hydraulic depth D and hydraulic radius R both vary 
h y. Thus, for Tolkmitt’s special case, the hydraulic exponents20 can be 
wnas N= M=4. 

For channel sections of a general shape, the flow-profile equation20 is 


N 
Yn me ste Ye uo yN-M 
x =——|u- +{— —— du|+constant. .(8 
So J oo J few (8) 


ie ees 


sre, using author’s notation, yp = d, So = S, u = (d + z)/d and (d + h)/d for a 
sing of water level, and u = (d - z)/d and (d - h)/d for a lowering of water 
el, When Chezy’s formula is “ised, it can be shown that yg/yn = sC2/g 1/4, 
is, Eq. 8 can be written partially in author’s notation as 


N/4 N-M 
d us du sC2 uiu 
= a ——— + —— du | + constant. . .(9 


re ix a = F(u, N), a varied-flow function. For Tolkmitt’s special 
u 


fa 
e, Eq. 9 is reduced to 
d eCA\ et (s du 
= — - - — teins, cha ane 10 

x : E ( rs 13, ; “=, |. constan (10) 
i u 2 1 
4 dus = A, util a a 

Sy see ="F (iu "4) = 5 nq) +5 tan u..... (11) 


lay be noted that F (u, 4) = 7/4 when u becomes infinity. 
for a raising of water level, X =u =(d+z)/d=1+z/dandu>1. Thus, 


11 becomes 


F(u, 4)= Gin (1429) +5 tant a+) sect a (12) 


May, 1960, by R. D. Goodrich. 
Prof, of Hydr. Engrg. in charge, Univ. of Illinois, Urbana, Il, 
Open-Channel Hydraulics,” by Ven Te Chow, McGraw-Hill Book Co., Inc., New 


1959, 
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From Eq. 2 . 
f(X) = u - F(u, 4) + y 


In the varied-flow function table,29 for u > 1, a constant, 1/4, was subti 
ed from F(u, 4) to make the function value near to zero when the varia 
becomes very large. Let F (u, 4) represent the table value. Then, 


{(X) =a =F (u, 4) <. «246 eae 


This relation indicates that the Tolkmitt’s table canbe replaced by the gel 
varied-flow function table for flow-profile-computation purposes. Fol 
ample, with u = 1.060, the table value29 for F(u, 4) for N = 4.0 is 0.506 
Eq. 14, £{(X) = 1.060 - 0.506 = 0.554. The above relation also applies to | 
= (d+h)/d=1+h/d. 

Similarly, the varied-flow function for a lowering of water level, with; 
= (d - z)/d = 1 - z/d and u <1, can be shown as 


F(u, 4) = 5 In @Z-+a tan’? (1 - 3) ae 


From Eq. 4, 
£(x) = F (u,4) a eee ee 


Thus, Tolkmitt’s table value is identical with the value given by the va 
flow function table. This also applied to X = u = (d - h)/d=1 - h/d. 

Tolkmitt’s table is widely known in European literature. The author s 
be commended for introducing this table to the American literature. 
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HOOD INLET FOR CLOSED CONDUIT SPILLWAYS2 


Discussion by C. D. Smith 


. D. SMITH, 19 F. ASCE.—It is heartily agreed that the hydraulics of cul- 
3 is not as simple as was once believed. In this regard, Mr. Blaisdell’s 
r represents a further useful contribution to this interesting and important 
apt. 
ne feature of the hood inlet which was not mentioned, and which should re- 
> atiention, is a discussion of the limitations of the design. In the writer’s 
on, there are very definite structural and hydraulic limitations which 
ld be considered by the designer when he is making his selection of inlet 
mn. 
he first point concerns the design of the splitter-type anti-vortex wall. In 
, a design using the hood inlet was carried out under the direction of the 
r for a 60-in. corrugated iron pipe. At that time, the only data reported 
he anti-vortex device was for the splitter-type wall. According to the 
mum dimensions, the splitter wall would have to cantilever upstream 
) ft beyond the invert of the pipe. It was found that several ring stiffeners 
d be required for the pipe, and structural steel angles would be required 
inforce the thin splitter wall. The design was complicated, and there was 
some question as to the adequacy of the splitter wall to resist impact due 
dating ice cakes or debris. Finally the design was abandoned in favor of 
re conventional reinforced concrete inlet with a rounded lip. 
may be that problems similar to this were the reason for the objections 
d to the splitter type anti-vortex wall. The anti-vortex plate of Fig. 8(f) 
ated to be equally efficient. It appears to be much simpler structurally, 
hould be a satisfactory replacement for the wall. 
1@ second point concerns the limitations of any design which gives slug 
Slug flow is unsteady flow and associated with such flow are varying 
sures, noise, vibration, and pounding at the inlet, and intensified wave ac- 
it the outlet. It appears that acceptability or rejection of these flow con- 
1S depends on the size of the works. For small stock-watering dams us- 
n 18 in. or 24 in. pipe these problems are probably not critical. For large 
of 60 in. or greater it would seem a wise precaution to avoid designs 
1 give slug flow. 
Should be mentioned that when the spillway operates, slug flow would pre- 
nost of the time. This is because the inlet primes with a fairly low dis- 
e, but the pipe does not flow full continuously until this discharge has 
exceeded several times. 
@ final and most important point concerns the development of sub- 
spheric pressures in the inlet. It was stated under the heading “Theory” 


fay, 1960, by Fred W. Blaisdell. 
agnoc. Prof. of Civ, Engrg., Univ. of sisiidicten: Saskatoon, Canada, 
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that if the pipe is steep, the pressure “will ordinarily be sub-atmosphi 
throughout most of the conduit length.” This would also be true for any 0 
type of inlet on a steep pipe flowing full. However, the local contractio 
flow in the hood inlet will give a pressure much lower at that point than é 
where in the pipe. 

Consider a sharp edged re-entrant tube flowing full: The coefficient of 
jet contraction C¢ = 0.5, and the coefficient of head loss 


2 
Ke “(es S ) = 1.0. 


The pressure head at the invert at the inlet will be 
P v2 


Ww 2¢g 


where V is the local velocity of the contracted jet, and not the average vel 
in the pipe. Therefore, 


J. te His 
Ww 


With C, = 0.5, the local velocity head will be 4 times as great as the ve: 
ty head based on average pipe velocity, and the pressure head will be co 
spondingly lower. Eq. 16 has been verified experimentally by the write 
others for a sharp edged re-entrant tube flowing full. 

Since the value of Kg for the hood inlet is also close to unity, one migh 
pect that C, would also equal 0.5. Actually, the cases are not:exactly ¢ 
parable, as the bevel cut on the hood inlet has some modifying influence. 
not possible to determine Cg for the hood inlet flowing as an orifice, bec 
it will not flowas an orifice. However, the “effective” value of Cc, can be 
culated from Eq. 16 if Q, A, H, and P/w are measured. . 

In tests conducted by the writer, piezometer readings were taken alon 
pipe invert downstream from the hood inlet. The pressures were not a 
for the hood inlet as for the sharp edged re-entrant. It was calculated thé 
effective value of Ce was approximately 0.6. This is a higher coefficient 
for the ordinary sharp-edged re-entrant, but it appears reasonable fo: 
case. 

It was also observed that the crown pressures along the hood inlet 
greater than 1 diameter above the invert pressures. This also appears 
sonable when the upward component of velocity, directed toward the cro 
the inlet, is considered. Actually, the flow pattern is exceedingly con 
and the value of Cc = 0.6 is suggested only as a value for computing press 
rather than giving an actual area of flow in the inlet. 

Eq. 16 may be written in dimensionless terms by dividing by D 


ean 
Pw H _\CcA 


D D 2¢D Se) Leite #, Bie abe Pel ge) epee ee * 


Substitution of Ce = 0.6, and reduction, yields an equation inthe same’ 
as Fig. 6: : 
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) illustrate the significance of Eq. 18, suppose a prototype pipe 48 in. in 
ster were to be built with such a length, slope, and roughness as to yield 
ischarge curve shown in Fig. 6 for a pipe slope of 0.200. This pipe should 
ull when H/D = 2 and give a discharge corresponding to Q/D9/2 = 15.5. 
itution of these values in Eq. 18 shows that the pressure head at the invert 
{be -14.8 pipe diameters. This would correspond to -59.2 ft on the pro- 
2, SO obviously it would not perform as indicated onthe model under these 
tions. 

yalue of -1/2an atmosphere is often imposedas adesign limit in hydrau- 
ructures. Lower pressures are considered to be in the realm of potential 
ution. If this value was used in Eq. 18, then the permissible maximum 
9f pipe for the conditions of the above example would be 14 in. 

the writer’s view, the negative pressures at the inletimpose a very defi- 
mitonthe size of pipe for which the hood inlet may be used with success. 
jesigner should check each particular case, using Eq. 18, to insure that 
ures are within acceptable limits. 
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p- 95. Delete lines 5 and 6 of paragraph 3 and substitute therefore: 
had proved, the dissipation of energy balances nearly with the productior 
in the range of approximately 0 < y/h < 1/2, and with the diffusion rate - 
outer region (0.8 < y/h < 1), and the turbulence 


p. 95. Delete Eqs. 25 and 26 and replace with 


eto 
Dissipation = A’ pv ae A p\|———] ~~ se ee 
N Ue - 
ee ee: waeaw) bit u2 vy 
Diffusion = ay (o u 5} erat | rag pP =} 
Sa pu,” d u2 v 
1 h d(y, h U3 ee @©688 © 6 @€ @ @ Cie 


Considering the isotropic character of turbulence in the outer (open chi 
or center (pipe) region, Eq. 25' can be written as 


pU,3 
Fras hares sae eee ; 
Therefore, by equating Eq. 25 to Eq. 24 and Eq. 25'' we have 


Diffusion = Ay 


3 
Jn (Ma2\ _ 
Sea = £1 (y/h) 


(0 < y/h < 0.5) 


ae 3 as testis bi Gate eae a 
B' |= = f , (y/h) 
U, | 


(0.8 < y/h< 1) 
p. 96. Just prior to the last paragraph insert 


“= 


JO a 
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», the diffusion coefficient for uniform open channel flows is givenfrom 
) as 
Fe Dz, = const. Vg h3/2 gel/2. (31) 
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In Eq. 9, insert a parenthesis before the term 0.15 and after the term 
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